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Various works performed by the present authors in the 1990s are reviewed. The topics 
discussed in this paper are mainly related to the time-evolution of the coherent and the 



, squeezed states of the systems obeying the su(2)- and the s?i(l, l)-algebra. The formulations 



are based on the basic idea of the boson mapping and the TDHF theory in canonical form. 
Under the time-dependent variational procedure for trial states appropriately chosen, the 
time-evolution of the systems under investigation is described. Further, it is shown that this 
method enables us to obtain the classical counterparts of the original quantal systems. 

§1. Introduction and preliminaries 



o 

00 , 

. For the microscopic studies of nuclear collective motions^ the year 1960 is unfor- 

gettable. In this year, a theory was proposed by Marumori,tP also independently by 



Arvieu &: Veneroni and Bar anger. © This theory is called the quasi-particle random 
phase approximation (QRPA) and, with the help of this theory, microscopic struc- 
tures of collective vibrational states, especially, the first excited states of spherical 
^ ■ even nuclei, were well described. Under the success of the QRPA theory for the first 

. excited states, the next was a problem how to describe the higher excited states, for 

example, the well-known (0^, 2+, 4"*") triplet observed in the second excited states. In 
vj ■ response to the above-mentioned situation, in 1962, Belyaev & Zelevinsky proposed 

. so-called boson expansion theory© and slightly later, in 1964, Marumori proposed 

I also boson expansion theory, together with Yamamura (one of the present authors) 

and Tokunaga.o' The former and the latter are called the BZ and the MYT method. 
Both are suitable for the case of fermion-pairs. The case of particle-hole pairs was 
formulated by da Providencia (one of the present authors) &: Weneser and Mar- 
shalek.EP'EF These are called the boson expansion theory all together. In Ref.^), the 
relation between the BZ and the MYT expansion was investigated from the side of 
the former form. We can find various further studies concerning the boson expansion 
in the review by Klein & Marshalek in 1991. & 

However, we should draw a distinction between the BZ and the MYT method. 
The former is based on the boson realization of fermion-pair algebra (Lie algebra) 
and the latter based on the mapping from fermion- to boson-space, which is called 
the boson mapping. The boson mapping is independent from the Lie algebra and. 
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if the Lie algebra is treated in the framework of the boson mapping, we can arrive 
at its boson reahzation. In this sense, it should be noted that the boson mapping is 
conceptually wider than the boson realization of the Lie algebra. However, it is not 
necessary for the boson mapping to restrict ourselves to the fermion- and the boson- 
space. For example, let us consider a quantal system which can be described in the 
framework of the orthogonal set {|n));n = 0, 1,2, ■■■ ,A''}. It is not always many- 
fermion system. Next, we take up another set which is not always composed 

of the boson operators, but, here, following the boson mapping, composed of a kind 
of boson operator (a, a*) : 

In) = 4=r(«*riO) • (1-1) 



'77,! 

Then, we set up the one-to-one correspondence between the two sets : 

|n))~|n). (n = 0,1, 2, •••,7V) (L2) 
In the boson mapping, the following operator Um plays a central role : 

N 
n=0 

With the use of Un-, we have 

UN\n)) = \n), (n = 0,1, 2, •••,7V) (L4) 

uUr^) = \n)), (n = 0,1, 2, •••,7V) (L5) 

ul,\N + k) = Q. {k = l,2,Z,---) (L6) 

Usually, the space spanned by the set {|n); n = 0, 1, 2, • • • , TV} is called the physical 
space. The operators JJn and tf]^ obey 



u],Un = 1 , (1-7) 

N 

UnUI = Pn , Pn = Y1 I") ("I ■ (1-8) 



n=0 



Here, Pn is a projection operator to the physical space which satisfies Pj^ = Pn and 
P^ = Pn- Next, we consider the operator Oi, in the original space which satisfies 



{{n + iy\0,\n)) = J ———0,{n) . (z/ > 0, n = 0, 1, 2, • • • , - i.) (L9) 
y (n -I- uy. 

Here, Oy{n) is a function of n which is obtained by calculating the matrix element 
((n -I- v\Oi,\n)). Then, in the boson space, O^, is transformed to UnOuUn ^^'^ 
expressed as 

UnO.uI = PnO.Pn , (MO) 
= {aYOuin) , h = a*a . (Ml) 
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The above is an outhne of the boson mapping and in later section, we will use the 
above-mentioned scheme. 

On the other hand, there exists another powerful approach to the microscopic 
theory of collective motion, which has been called the time-dependent Hartree-Fock 
(TDHF) theory. At the early stage, the idea of the TDHF theory was proposed 
by Nogami® and slightly later by MarumoriEl^ for the case of the small amplitude 
vibrational motion. The frequency of the small fluctuation around a static HF field 
is calculated in the frames of the TDHF theory. At the middle of the 1970s, the 
TDHF theory was revived not only for the study of anharmonic vibration but also 
for the studies of heavy- ion reaction and nuclear fissions. Many ideas or forms were 
proposed for various problems by many authors. lP Especially, in 1980, a canonical 
form for the TDHF theory was proposed by Mammori, together with Maskawa, 
Sakata and Kuriyama (one of the present authors). £3-* The TDHF theory is a possible 
method for obtaining an approximate solution of the time-dependent Schrodinger 
equation. First, as a trial state for the time-dependent variation, we prepare a Slater 
determinant under a definite form containing parameters. Through the variational 
procedure, we obtain certain differential equations with respect to t and, then, by 
solving under an appropriate initial condition, the time-dependence of the parameters 
is determined. Thus, we have an approximate solution of the original time-dependent 
Schrodinger equation. In the canonical form of the TDHF theory, the parameters 
are required to obey certain relation which the present authors call the canonicity 
condition. Then, the parameters can be regarded as canonical variables in classical 
mechanics and the certain differential equation with respect to t is reduced to the 
Hamilton equation of motion in classical mechanics, the Hamiltonian of which is given 
as the expectation value of the starting quantal Hamiltonian for the trial function. 

Introduction of the canonicity condition into the TDHF theory suggests us to 
formulate the variation for the time-dependent Schrodinger equation in relation to 
classical mechanics for any quantal system. Let H denote the Hamiltonian under 
investigation, which is not always many-fermion system. For H, we set up the time- 
dependent Schrodinger equation 

d 

z^— |exact) = ff|exact) . (1'12) 

As for an approximate solution of Eq. ( |1-12| ), we prepare a trial state \p,q). Here, p 
and q denote two real parameters obeying the relation 

{P,q\in'—\p,q} =P+-Q^ ^ 

{p,q\ih-^\p,q) = ^ . (1-13) 
op op 

Here, S is arbitrary real function of p and q. The present authors call the relation 
(1-13) the canonicity condition. Then, we have the relation 

d dS 
{p, qlif^Q^ - H\p, q) =pq- H{p, + > (1-14) 

H{p,q) = {p,q\H\p,q) . (1-15) 
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The relations (1-14) and (1-15) can be regarded as the Lagrangian and Hamiltonian, 
respectively and the variation gives us 



j {p, q\in^^ - H\p, q)dt = sj{pq- H{p, q))dt = . (1-16) 



Thus, we have 



— p--— (117) 
dp ' dq 



The differential equation is nothing but the Hamilton equation of motion in classical 
mechanics. Therefore, (p, q) can be regarded as canonical variables. By solving Eq. 
( 1-17| ) under an appropriate initial condition, p and q can be determined as function 



of t. The substitution of the solution of Eq. ( |1-17| ) into \p, q) gives us an approximate 
solution of Eq. (1-12). If S is chosen as the form S = pq/2, the canonicity condition 





is written as 

X* X 

{p,q\dx\p,q) = — , {p,q\dx*\p^Q) = -y > (1-18) 

(1-19) 

Of course, {X,X*) denote the complex parameters. The above is an outline of the 
variational method for the time-dependent Schrodinger equation which is suggested 
from the canonical form of the TDHF theory. In later section, we use the above- 
mentioned scheme. 

The above-mentioned time-dependent variational approach to quantal system 
has two characteristic aspects. One is to give an approximate description of time- 
evolution of the quantum state. The other is to give a possible classical counterpart 
for an original quantal system under a suitable choice of the trial state. Then, it is ex- 
pected that the original system is reproduced in a disguised form under an appropri- 
ate requantization procedure. This viewpoint was originally proposed by Marshalek 
and Holzwarth for the case of the time-dependent Hartree-Bogoliubov theory. EIP 
They obtained the Holstein-Primakoff boson representation for the fermion-pair al- 
gebra, which is of the same form as that given by the boson mapping. B> 

With the aim of formulating the TDHF theory and its extension in the canonical 
form, the present authors (A. K. &; M. Y.) have investigated the above-mentioned 
idea extensively for the cases of various many-fermion systems. For example, with the 
aid of the Grassmann variables, odd fermion system can be treated in the framework 
of the above-mentioned scheme. These works have been reviewed in the paper by 
the present authors (A. K. & M. Y.). Further, the above idea was applied to 
the investigation of the behavior of the systems in highly excited states also by the 
present authors (A. K., J. da P. & M. Y.).ll3) The basic idea for this case is in 
the description of many-fermion systems under the phase space doubling, which is 
fundamental in the thermo field dynamics formalism. US This description gave us 
the so-called thermal boson expansion and the time-evolution of the thermal states. 
The above all papers are related with many-fermion systems. 



Time- evolution of the coherent and the squeezed states 



5 



Focussing on the su{2)- and su{l, l)-algebra, in this paper, we will review some 
investigations performed mainly by the present authors after 1991. Of course, the 
topics discussed in this paper are related to the title of this paper. It is well known 
that the su(2)-algebra has played an important and interesting role in the study 
of many-body physics. Even if restricted to nuclear theory, this algebra has served 
us not only for schematic understandings of collective dynamics such as the Lipkin 
model but also for studies of realistic phenomena such as superconducting phase and 
rotational motion. We can learn various investigations related to the ■SM(2)-algebra 
until 1990 in Refs.|7[) and 12). A possible boson representation of the SM(2)-algebra 
was given by Schwinger and this form will play a central role in this paper. 
On the other hand, as was originally discussed by Celeghini, Rasetti, Tarlini and 
Vitiello,!!^ the su{l, l)-algebra is alsomiite interesting. With the aid of this algebra 
in the Schwinger boson representation, Ej) we can describe the "damped and amplified 
oscillator" quantum mechanically in the conservative form. The investigation of the 
behavior of systems in highly excited states, such as nuclear phenomena observed 
in highly excited states, is one of the interesting problems in quantum many-body 
physics. In the theoretical treatment, it is necessary to deal with the systems at 
finite temperature. The above means that the su{l, l)-algebra presents us a possible 
entrance to the problems related to finite temperature. The interest of the present 
authors in the 1990s has been concerned mainly with the time-evolution of many- 
body systems obeying the su{2)- and su{l, l)-algebra. Of course, as is clear from the 
title of this paper, the interest is related with the classical descriptions of quantum 
systems, and the technique of the boson mapping and the time-dependent variational 
method are fully used. These investigations are divided into two groups. One is 
ated mainly with time-evolution of the coherent states in many-boaon svstems 
a) and the other with the squeezed states in jnamir-body systems, eIP "E^ some 
of which were published together with Fujiwara. E^' e2P Our interest is concerned 
not only with the coherent and the squeezed state for a simple system but also with 
those for the systems obeying the su{2)- and the su{l, l)-algebra. Starting from the 
conventional boson coherent state, we formulate various forms which are suitable 
for treating the su{2)- and the su{l, l)-boson models. Especially, one of the present 
authors (Y. T.) investigated the formalism, which is partially presented inthis paper, 
in relation to the WKB result and the Maslov phase in the limit ofh — > 0. E3) Modified 
forms of the conventional Holstein-Primakoff boson representation of the su{l, 1)- 
algebra and g-boson realizations of the sUq{2)- and the suq(l, l)-algebra, which we 
will not contact with in this paper, are discussed in Refs. ^ and ^), respectively. 

In §§2 and 3, starting from the Schwinger boson representation for the ■sm(2)- 
and the su{l, l)-algebra, the Holstein-Primakoff boson representations are derived 
under the use of the basic idea of the boson mapping. Further, with the aid of the 
coherent state, the classical counterparts of the both algebras are obtained. Section 
4 and 5 are devoted to investigating the "damped and amplified oscillation" and 
the sn(2)-algebraic model in the framework of the su{l, l)-algebra and its coherent 
state. The idea is found in the phase space doubling appearing in the thermo field 
dynamics formalism. In particular, the thermal effects observed in their systems 
are discussed. In §§6 ~ 9, a boson system interacting with the external harmonic 
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oscillator is described. Under a certain condition for the interaction, the results 
are almost the same as those obtained in §§4 and 5. In §10, the time-dependent 
variational approach with the squeezed state is presented and, in §§11 and 12, the 
squeezed state for the su(2)-algebraic models of both many-boson and many-fermion 
systems are discussed. Finally, the concluding remarks are mentioned in §13. 

§2. The sii(2)-algebra — The boson realization and its classical 

counterpart 

First, we will recapitulate the boson representation of the sii(2)-spin system in 
a form suitable for later discussion. For the sake of completeness, various familiar 
relations will be listed. Let us introduce a boson space constructed in terms of two 
kinds of boson operators (a, a*) and {b,b*). In this space, we define the following 
three operators : 

S+ = ha*b , S- = nb*a , = |(a*a - b*b) . (2-1) 

Here, h denotes the Planck constant. The above three operators S±fi form the 
sn(2)-algebra, which obeys 

5* = 5+ , = So , (2-2) 

[S+, S^] = +2hSo , [So, S±] = ±nS± . (2-3) 



The form (24) is called the Schwinger boson representation of the sti(2)-algebra.li^ 
The Casimir operator rgu{2) can be written as 



A.(2) = Si + \{S-S+ + S+S-) = S{S + h) . (2-4) 



2 

Here, S and its property are given as 

S = +'^{a*a + b*b) , (2-5) 

[5,S±,o]=0. (2-6) 

We can understand that S is an operator denoting the magnitude of the su(2)-spin. 

The eigenstate of S and 5*0 with the eigenvalues S and Sq, respectively, is ob- 
tained in the form 

\S,So) = ^ I,, (a*)-+-°(fe*)--°|0) , (2-7) 

v/(s + so)!(s - So)! 

a|0) = 6|0) = . (2-8) 

The eigenvalues S and 5*0 are expressed as 

S = hs , So = nso , (2-9) 
1 3 

•s = 0, -,1, -,• • • , So = -s,-s + 1,- • • ,s - l,s . (2-10) 
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We call the space spanned by the set {I^, S'o)} for fixed value of S the s'u(2)-space 
with the magnitude S and at several occasions, we will use the notation |s,,so) for 
\S,So). 

Following Refs.p3D ~ let us derive the Holstein-Primakoff representation c2P 
from the Schwinger representation. Basic idea comes from the MYT boson mapping 
method proposed by Marumori, Yamamura and Tokunaga. For the above-mentioned 
aim, we prepare a new boson space spanned by the set {|n) = [\/nI]~^ • (X*)"|0); n = 
0, 1, 2, • • •}. Here, {X, X*) denotes boson operator and |0) is the vacuum for (X, X*). 
First, following the MYT boson mapping, we set up the following correspondence 
between the original and the new boson space : 

|s,so) ~ |s + So) • (S = hs , So = hso) (2-11) 

The spaces spanned by the set {|s + sq); sq = —s, — s + 1, • • • , s — 1, s} and the set 
{|2s + k);k = 1, 2, 3, • • •} are physical and unphysical space, respectively. Under the 
correspondence ( 2T1[ ), the following mapping operator U is introduced : 



U= \s + so){s,so\ . (2-12) 

so=-s 

The properties of U are shown in the form 

s 

WU = \s,so){s,so\ = 1 , 

so=-s 
s 

UU^ = J2 \s + so){s + so\ = P2s ■ (2-13) 

so=-s 

Here, P2s denotes a projection operator to the physical space : 

^2. = ^2., Pl=P2s. (2-14) 

From the above conditions, the following relations are derived : 

U\s,so) = \s + so) , (2-15) 
U'^ls + So) = \s,so) , 

U''\2s + k) = . (/c = 1,2,3, •••) (2-16) 
With the use of U, the operators S±fi are transformed to the form 

US±nU^ = P2s • S±n ■ P2s , 



s+ = Vhx* ■ J2S-nx*x , 



s_ = ^2S -hx*x -Vhx , 

^0 = nX*X - S . (2-17) 
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The above operators S±fl satisfy 

Q C** Q 

0_ — , Oq — OQ , 

[s+, = +2nSo , [So, s±] = ±hs± 

So + \{S-S+ + S+S-) = S{S + h) . 



(2-18) 
(2-19) 

(2-20) 



The relations (2T8) ~ (2-20) correspond to Eqs. (2-2) ~ ( |2-4| ). The above is nothing 
but a refined form given by MarshalekcP and it is called the Holstein-Primakoff boson 
representation. 

Ou]i_thini problem is concerned with the coherent state I7) for the ■su(2)-spin 
system. ESP For the construction of the state I7), we note the form 



\s,so) 



(2s)! V (s + so)! V ^ 
Then, we set up the following correspondence : 

1 



is-so)l (S^^ .(r)2^|0). 



(2-21) 



{s - So)! / S- 



S + So 



V S- 



(2s)! V (s + so)! 



— ~ exp — • — — = exp — • a*b 



u h 



{b*y'\0) - exp ( \ljj^wub*\ |0) . 



Here, w and v denote complex parameters and u is given as 



Under the correspondence (2-22), we define the state I7) in the form 



I7) = Mc • exp 



V Sj, 



■ exp 



u h 

Here, Mc denotes the normalization constant given by 

|2\ 



wuh* |0) . 



Mc = exp 



\w\ 



(2-22) 



(2-23) 



(2-24) 



(2-25) 



We can see that in the state I7) the part exp ( \/2/ hwub* ) 1 0) is a superposition of the 
states with the maximum weight (the operation of S- on this part makes it vanish) 
and the part exp[(f/u) ■ S+/h] denotes a superposition of the raising operation. The 
state 1 7) can be rewritten in the form 



It) = exp 



\w\ 



■ exp 



—w{va* + ub*] 



|0) 



It may be clear that I7) is a boson coherent state, which satisfies 



a|7) 



wvl'j) , 6I7) 



(2-26) 



(2-27) 
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In other representation, we have 

{ua - vb)\j) = . (2-28) 

Here, v*a + ub and ud — vb represent annihilations of independent boson operators. 
The expectation values of S and 5-1-^0 ^-re calculated in the following form : 

(7|^|7) = {S)c = \w\^ , (2-29) 
(7|S+|7) = (5+), = 2|«;|Vu, 
(7|S_|7) = = 2\w\^uv , 

{^\Soh) = {So)c = \w\'{\v\'-u') . (2-30) 

Instead of the parameters {w,w*) and {v,v*), we introduce the parameters ("0,5^) 
and {X,X*), which satisfy the conditions 

i^lihd^h) = s , {j\inds\^)=o, (2-31) 

(7|5xl7) = +4^, {l\dx*h) = -§- (2-32) 

We can learn in the TDHF theory that the parameters.^ (■0, S) and {X,X*) are 
regarded as canonical variables in classical mechanics. liS^'EJ) The variables {tp, S) are 
called the phase angle and the action variable and (X,X*) are classical variables in 
boson type. For S, we set up the relation 

(5), = = S . (2-33) 

The conditions (2-31) and (2-32) give us the following relations : 

w = \/S exp f"^"^ 




,^X, u=^l-^X*X. (2-34) 

Substituting the relation (2-34) into the expectation values {S±fl)c in Eq. (2-30), we 
have 



= Vnx* ■ V2S - hX*X , 

= V2S - hx*x ■ Vnx , 

{So)c = hX*X - S . (2-35) 

The expectation values (5'±^o)c shown in the relation (2-35) satisfy 

= , (So): = {So)c , (2-36) 

[(5+)„(5_),]p = H)[2(5o)e] , 

[iSo)c,iS±)c]p = i-i)[±iS±)c] , (2-37) 

{So)l + liiS-US+)c + iS+US^)e) = . (2-38) 
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Here, [A, B]-p denotes the Poisson bracket defined by 

fdAdB dBdA\ 1 / dA dB dB dA . 



dip dS dipdSj ih \dX dX* dX dX' 

The relations (2-36) ~ (2-38) correspond to Eqs. (2-18) ~ (2-20) and, in this sense, 
the set {(5'±^o)c} can be called the classical counterpart of the set {S'i^o} and, further, 
{S'i^o}- Under the replacements {X,X*) {X,X*) and S ^ S, the classical 
counterpart becomes the quantal form. 

§3. The su{l, l)-algebra — The boson reahzation and its classical 

counterpart 

Main aim of this section is to recapitulate the boson representation of the 
su{l, l)-spin system in a form suitable for later discussion. In this case, various 
relations are not so familiar to those in the S'u(2)-spin system. In a manner similar 
to the case of the ■sn(2)-spin, we define the following three operators : 

f+ = ha*b* , f_ = hba , Tq = ^{a*a + bb*) . (3-1) 

The operators T±fi form the su{l, l)-algebra, which obeys 

fl = f+, To* = To , (3-2) 
[f+, f_] = -2hfo , [To, f±] = ±nf± . (3-3) 

The form ( |3T| ) is called the Schwinger boson representation of the su{l, l)-algebra. 
The Casimir operator /^^^(i^i) can be written as 

rsu(i,i) = To - liT-f+ + r+r_) = f (f - h) . (3-4) 

Here, T and its property are given as 

f = -^{a*a-bb*) , (3-5) 

[f,T±,o]=0. (3-6) 

The operator T is not positive-definite, but, in a sense analogous to the case of the 
■SM(2)-spin, T may be called the magnitude of the su{l, l)-spin. 

The operator T is not positive-definite and the eigenvalue of Tq is larger than h/2. 
Then, compared with the case of the su(2)-spin system in the relations (2-7) ~ (2T0), 
the treatment of the eigenstate of T and Tq is rather complicated. The eigenstate 
of T and Tq with the eigenvalues T and Tq, respectively, is formally obtained in the 
form 

|r,To) = ==(r)^o'*(b-)*°+*-^|0) , (3-7) 

T = nt , To = fito . (3-8) 
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Since to > t, tQ > —t + 1 and to ^ 1/2, the state \T,Tq) is divided into two cases : 

(i) t = i 1,|2,---, to = t,t + l,t + 2,--- , (3-9) 

(ii) t = 0, -1, • • • , to = -t + l,-t + 2,-t + 3,--- . (3-10) 

In this paper, we restrict ourselves to case (i), i.e., the case t > 1/2. We call the 
space spanned by the set {\T,Tq)} for fixed value of T (> fi/2) the su(l, l)-space 
with the magnitude T and at several occasions, we will use the notation |t,to) for 
\T,To). 

In the same manner as that in the case of the s'u(2)-spin, let us derive the 
Holstein-Primakoff representation. IIZP We set up the following correspondence be- 
tween the original and the new boson space : 

\t,to) ^\to-t) . {T = ht , To = nto) (3-11) 

Since t^ = t,t + l,t + 2, ■ ■ the whole space in the new space is physical. Under the 
correspondence ( 3-11 ), the following mapping operator V is introduced : 

E l*o-i)(i,io| • (3-12) 

to=t 

The properties of V are shown in the form 

V^V=Y,\t,to){t,to\ = l , 

to=t 

VV^ =Y,\to-t){to-t\ = l . (3-13) 

to=t 

In contrast with the relation (2-13), there does not exist the projection operator such 
as in the relation (3-13). From the above conditions, the following relations are 
derived : 

V\t,to) = \to-t) , (3-14) 
V^\to-t) = \t,to) . (3-15) 
With the use of V, the operators T±fi are transformed to the form 



r+ = Vhx* ■J2T + hx*x , 



r_ = V 2T + hx*x ■ Vhx , 

fo = hX*X + T. (3-16) 
The above operators T±fi satisfy 

f * = f+ , To* = To , (3-17) 
[r+, f_] = -2nfo , [To, T±] = ±hf± , (3-18) 
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The relations (3-17) ~ (3T9) correspond to Eqs. (3-2) ~ ( p-4D . Judging from the 
form, it can be called the Holstein-Primakoff boson representation. 

In parallel with the case of the sti(2)-spin system, our third problem is concerned 
with the coherent state |c) for the 5^(1, l)-spin system. HIP For the construction of 
|c), we note the form 

|t,to) = ^= ^{'^ •(r)2*-^|0). (3-20) 

Then, we set up the following correspondence : 

1 ff.V"'' (V fA (Y 



— ~ exp — • — = exp \ — -ah 



V(*o-t)!(*o + t-l)! V ^ / ^/ Vf/ 

(6*f-i|0)~exp(^y|^6*^|0) . (3-21) 

Here, W and Y denote complex parameters and U is given as 



\J=^\^\Y\^. (3-22) 
Under the correspondence (3-21), we define the state |c) in the form 

|„=«,.expg.|).exp(;ffr)|0). (3.23) 

Here, denotes the normalization constant given by 

1 / \W?\ 

In this case, also, we can see that in the state |c) the part eyj^{\plJfiW jU • 6*)|0) is 
a superposition of the states with the maximum weight (the operation of r_ on this 
part makes it vanish) and the part exp[(l//[/) • T^jfi\ denotes a superposition of the 
raising operation. The state |c) satisfies 

{Vh-Ya*)\c) = ^^\c) , 

(JJa-Y}f)\c) =^ . (3-25) 

The operators Uh — Ya* and Ua — Yh* represent annihilation of independent boson 
operators and |c) is a coherent state in the sense given in the relation (3-25). As 
is clear from the form ( 2-26| ), the state I7) is one- mode coherent state and, in the 
present case, |c) is mixed- mode coherent state. 

The expectation values of T and T±i^ are calculated in the following form : 

(c|f|c) = (r), = |M^|2 + ^, (3-26) 
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(c|f+|c) = (r+), = 2(|Vr|2 + |) v*u, 

2, 

{c\To\c) = (ro)e =[\W\-' + ^) + . (3-27) 



{c\T_\c) = iT.), = 2{\W\^ + ^) UV 



In the same manner as in the case of the su(2)-spin system, we introduce the pa- 
rameters (0,T) and {X,X*) for the old {W,W*) and {V,V*). For these variables, 
we set up the conditions 

{c\ind^\c) = T - ^ , {c\indrin\c) = , (3-28) 
{c\dx\c) = +^. (c|9x*|c> = -f . (3-29) 

The set of the relations (3-28) and (3-29) is called the canonicity condition. 
For T, we use the relation 

(T), = |T^|2 + ^ = T. (3-30) 
The conditions (3-28) and (3-29) give us the form 



Substituting the relation (3-31) into the expectation values (I±^o)c in Eq. (3-27), we 
obtain 

(r+)c = \/hX* ■ y/2T + hx*x , 
(r_)c = y/2T + hx*x ■ Vhx , 

{To)c = hX*X + T . (3-32) 
The expectation values (7±_o)c shown in the relation (3-32) satisfy 

(r_): = (T+), , (To): = (To)c , (3-33) 

[(r+)„(r_yp = H)[-2(ro),] , 

[{n)c,{T±)c]p = {-i)[±{T±)c] , (3-34) 
infc - \{{T^)c{T+)c + {T+)c{T-)c) = . (3-35) 

Here, [A, B]p is the same as that shown in Eq. ( p-39 ) under the replacement {tp ^ cj> 
and S — > T). The relations (3-33) ~ (3-35) correspond to Eqs. (3-17) ~ (3T9) and, 
in this sense, the set {(T±^o)c} is the classical counterpart of the set {T±fi}. This is 
in the same situation as the case of the SM(2)-spin. 
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§4. The sn(l, l)-algebra — One-dimensional "damped and amplified 
oscillator" and the sM(2)-algebraic model 

Following Ref.^), let us consider two-dimensional space, in which we define the 
following Lagrangian : 

L = mXY - kXY + m-f{XY - XY - 'yXY) . (4-1) 

Here, {X, Y) denote two kinds of real variables and m, k and 7 are positive constants. 
The Lagrange equations of motion lead us to the equations of motion 

mX = -{k + m7^)X - 2m-fX , 
mY = —{k + ?Ti7^)y + 2m'yY . 

The solutions of Eqs. (4-2) and (4-3) are expressed in the form 

X = exp(-7t) cos{ojt + a°) , 
Y = B^ exp(+7t) cos(a;t + , 

k 
m 



(4-2) 
(4-3) 

(4-4) 
(4-5) 

(4-6) 



The solutions (4-4) and (4-5) show the "damped and amplified oscillation." The 
Lagrangian (|4-1|) can be rewritten as 



L 



1 



-mx 



-kx' 



2^ -2^y 



m 



j{xy - xy) + ^(x^ - y^) 



Here, x and y are defined as 

1 



(X + Y) 



(X-Y). 



V2' " " ^/2 

The Lagrangian ( [4-7[ ) gives us the following Hamiltonian : 



1 



TT ( Px , ^ 2 2 

H = — muj X 

\ 2m 2 



^ + ^mcj^y^ ) - -f{xpy + ypr,) 
2m 2 I 



dL 

Px = -TT— = mx + m^y 
ox 



Py 



dL 
dy 



-my — m^x 



(4-7) 



(4-8) 



(4-9) 



(4-10) 



After canonical quantization, i.e., x ^ x, y ^ y, px ^ Px and py ^ py : [x,px] - 
[y,Py] = ih and [ other combinations ] =0, we obtain the quantized Hamiltonian : 



H = fiuob*b — fLOja*a — ihj{a*b* — ha) 



2uj T 



(4-11) 



Here, (a, a*) and {b,b*) are defined as 
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We can see that the Hamiltonian (4-11) is expressed in terms of T and T±. This 
means that the quantized "damped and amplified oscillator" can be described in 
the framework of the su{l, l)-algebra, in other word, it is a kind of the su{l, l)-spin 
system. 

The above treatment is interpreted as follows : In order to treat the system 
such as the "damped and amplified oscillator" as an isolated system, so-called phase 
space doubling is required. Then, the original intrinsic oscillator expressed in terms 
of the boson (6, b*) and the "external environment" expressed in terms of the boson 
(a, a*) appear and the interaction between both system is naturally introduced This 
is our interpretation, which permits us to generalize the Hamiltonian (4-11). Let 
the Hamiltonian of the intrinsic system be K given in the form 

K = F{nb*b) . (4-13) 

Here, F is a function of hb*b. Then, a possible generalization is performed by setting 
the following Hamiltonian : 

H = F{hb*b) - F{ha*a) - ihj{a*b* - ba) . (4-14) 

If F{hb*b) = hujb*b, the Hamiltonian ( [4-14| ) is reduced to the Hamiltonian (4-11). 
The Hamiltonian can be expressed in terms of the operators (T, T±o) : 

H = F{fo + f-h)- F{fo - f ) - ij{f+ - f_) . (4-15) 

For the convenience of late discussion, we rewrite the Hamiltonian (|4-15| ) in the form 

H = F{f,+f-fi)-F{f,-f) + 2-ify. (4-16) 

Here, we adopt the following form : 

f, = ^(f++f_), f, = -^(f+-f_) , t = fo. (4-17) 

We can find two ■SM(2)-algebraic models in nuclear physics. One is a model which 
enables us to describe pairing correlation in single-orbit shell model and another is 
for describing particle-hole correlation in two single-particle levels. In the first case, 
the intrinsic Hamiltonian K is given by 

k = 2e{So + S)-GS+S- . {S = n{2j + I) / A) (4-18) 

The quantities e and G denote the energy of the single-particle state and the strength 
of the pairing interaction, respectively. Further, j denotes the magnitude of the 
angular momentum of the orbit. Of course, the set {S±fi) forms the SM(2)-algebra 
and S denotes the magnitude of the sti(2)-spin. Then, if S±fl are replaced with S±fl 
shown in Eq. (2-17), K becomes of the form 



K 



G{S + \ 



{hb*b) + G-{nb*bf . (4-19) 
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Here, {X,X*) in Eq. (2-17) is replaced with {b,b*). Then, the Hamiltonian is given 
in the form 

H = A[e- G{S + h)] .(f-^^+4G-(f-^'jt + 2^fy . (4-20) 

The other is cahed the Lipkin model and the intrinsic Hamiltonian K is expressed 
as 

K = 2e{S + S)-^iSl + Sl) . (4-21) 

Here, S denotes the maximum magnitude of the SM(2)-spin that is equal to fifi/2, 
where fi is the degeneracy of the single-particle level. And e and G (> 0) represent 
the single-particle energy-spacing and the strength of the interaction, respectively. 
The set {S±fl) denotes the SM(2)-spin. In this case, together with the form (2-17), we 
adopt the basic assumption of the ATDHF theory and the semi-classical approxima- 
tion. Then, in the T-conserving framework, the Hamiltonian is given in the following 



H 



6hG 



e + G S + 



6G 



Here, ujq is defined by 



e + G S + 



T -- 
2 



2 U + G S + 



The quantity vq is a solution of the cubic equation 



4|6 + G|S+5 



n 

2 

1^0 



uo - 3hG = . 



(4-22) 
(4-23) 
(4-24) 



The detail of the above derivation has been given in Ref j20|) . The Hamiltonians 
( [4-20| ) and (4-22) depend on (T — fi/2), {T — h/2)Tz and Ty in the same form as each 
other and if the coefficient of the term (T — 7i/2)T^ vanishes, both Hamiltonians are 
reduced to the Hamiltonian (4T1). Hereafter, we will treat both systems in a unified 
way : 



H = 2ujT-[T 



+ 2/t 



T, + 27T, 



y ■ 



(4-25) 



The meanings oioJT and fx may be self-evident. Hereafter, we will omit the subscript 
"T" in OJT and /t. 

We are now at stage to discuss the time-evolution of the coherent state ( |3-23| ) 
for the Hamiltonian (|4-25D under the variational procedure. The expectation value 
of the Hamiltonian ( [4-25|) for the state ( |3-23|) , which we denote is expressed in 
the form 

h 



H = 2uj\T 
A 



+ 2\T^ + 2-^T, 



y ' 



A(T) = /-(1 + — 



T - - 

2 



(4-26) 
(4-27) 
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The expectation values of Tx, Ty and Tz for the state ( |3-23| ) are given as fohows : 

= ^{X* + X) ■ V2T + hX*X , 

Ty = -i^{x* - X) ■ ^2T + nx*x , 

Tz = hX*X + T . (4-28) 



Here, we used the definition (4-17) and the form (3-32) under the omission of the 
symbol ( )c- The set (4-28) satisfies the relations 

[Ty,Tz]p = Tx , [r„ Tjp = Ty , [Tx, Ty]p = -Tz , (4-29) 
[T, r,]p = [T, Ty]p = [T, Tz]p = , (4-30) 
- T^ - T^ = . (4-31) 

Since T and (j) are mutually canonical, the equations of motion are obtained as 

t = [T,H]p = 0, (4-32) 

= [0, H]p = 2uj + ^(2AT, + 2-fTy) . (4-33) 

Since T is a constant of motion, we set up 

T = I . (4-34) 

The quantity / is determined at the initial time. Instead of examining the Hamilton 
equations of motion for {X , X*), let us consider the equations of motion for T^, Ty 
and Tz : 

tx = [Tx,H]p = -2\Ty - 2^Tz , (4-35) 
ty = [Ty,H]p = +2XTx , (4-36) 
tz = [Tz,H]p = -2jTx , (4-37) 

From these equations, we have 

fx = -4(A2 - ^^)Tx . (A = A(T)) (4-38) 

Clearly, the solutions of Eq. (|4-38| ) are classified into four cases : 

(i) A2 > (ii) X^<-f\ (iii) A = -7 , (iv) A = +7 . (4-39) 

The results of the solutions of Eq. ( [4 -381 ) and the related points were given in Ref.|2l|) . 
In this paper, only the results are summarized. Including Eqs. (4-32) and (4-33), we 
adopt the canonical transformation from {cj),T; X , X*) to {6,I;XjJ)- 

Case (i) : A2 > 72 (Aq = ^/>^^) 



Tx = -\/J^-I^smx , 
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r, = Ay7^cosx-f J, 

AO Ao 

T = I , 

(f) = e + 2tan~^ 

H = 2u;(^I-^^+ 2X0 J , 

X = 2Xot + xo, e = 2ut + 9o 



XqI + a J + 7 V J' - ^ X 
XI + XoJ 2 



(4-40) 



1 dX . ,^ 
^^smx, (4-41) 



(4-42) 
(4-43) 



Case (ii) : A2 < 72 (70 = ^ff^) 



T^ = -^J^ + I''sin\ix , 
X_ 

70 



Ty = ^j2 + l2coshx + — J 

70 70 

r, = T^J^ + l^coshx- -J , 



70 
T = I , 



70 



= + 2tan~^ 



joI-XJ~WJ' + I\ .X 

— tanh — 

XI + JO J 2 



H = 2u[I--J+ 270 J , 
X = 270* + xo , = 2iot + Oo . 

Case (iii) : A = —7 



T = I 



(f) = e + 2tan"^ 



Jx 



I + J 



H = 2uj\^I - -j + 2XJ , 

X = 2Xt + xo , e = 2ut + 00 



(4-44) 



^dA . 



2 j^sinhx , (4-45) 



(4-46) 
(4-47) 



(4-48) 

(4-49) 

(4-50) 

(4-51) 



Case (iv) : A = +7 



xJ , 
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(4-52) 



(4-53) 



H = 2u;[I--]+2jJ 



(4-54) 

X = 2jt + xo, e = 2Ljt + eo. (4-55) 
We can see that J is also a constant of motion. The variables (X, X^) are given by 



Vhx 



Vnx* 



(4-56) 



Then, the time-dependence of various quantities is determined as functions of t for 
four cases. For example, \V\'^ given in Eq. (3-31) is expressed in the form 



\V\ 



nx*x 



n~i2 I n~i2 
^ X ^ y 



2T 2T(T + T^) 
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(4-57) 



Here, the relation (4-31) was used. It may be interesting to see that there exist four 
phases. 



§5. Description of thermal effects in terms of the coherent state for the 

su(l, l)-algebra 

In §4, we discussed the sii(l, l)-algebraic model, which permits us to describe 
influences of the "external environment" on the intrinsic system under investigation 
such as the harmonic oscillator. As the result of the interaction between two systems, 
the intrinsic system becomes the "damped and amplified oscillator." Therefore, it 
may be interesting to investigate such influences in the language of the thermal 
effects. We know that, in the thermo field dynamics formalism, EJ^ statistically mixed 
state can be described by the conventional Schrodinger equation with a modified 
Hamiltonian in a space with phase space doubling. In order to specify mixture of 
the pure states, the auxiliary variables are introduced. The formalism of the present 
sm(1, l)-algebraic model resembles the thermo field dynamics formalism. The boson 
operators (6,6*) describe the intrinsic system under investigation and (a, a*) are 
introduced for the auxiliary variables. The coherent state |c) shown in Eq. ( |3-23| ) 
contains both boson operators and, in the sense of the variational procedure, satisfies 
the Schrodinger equation derived from the intrinsic Hamiltonian ( 4-13| ) through the 
form (4-14). Therefore, |c) can he regarded as the state which describes the mixed 
state for the present system. ll2P'll2P'Ei|)'E3) 

First, let us expand partially the exponential function in the mixed-mode coher- 
ent state |c) in Eq. ( 3-23| ). Since the vacuum state |0), which satisfies the relations 
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a|0) = 6|0) = 0, is represented as |0) = |0)b ® |0)a (a|0)a = 6|0)fe = 0), the state |c) 
is expressed in the following form : 

oo 

\c) = Y.r{n-c)\\n{c))®\n)) . (5-1) 

n=0 

Here, ||n(c)) and \n)) denote normalized states for b- and a-boson, respectively. They 
are defined together with r{n;c) as follows : 

||n(c)) =e*"'^/2iV(n;|<5|2)-iexp f-^^ -1=(6*)" exp(<56*)|0)b , (5-2) 



\n)) = -^{a*r\^)a , (5-3) 
r{n; c) = |y|";7-("+i) exp N{n; . (5-4) 

Here, V is defined as F = \V\ exp(zcr/2) and N{n] \S\'^) and S are given by 



A'(";l*l^) = E^(m' 

r=0 ^ ' 



Ln(— I^P) : Laguerre polynomial , (5-5) 

nil ' 



(5-6) 



We should note that the set of the states {|ra))} in Eq. (5-3) is an orthonormal one 
satisfying ((m|n)) = 6mn- However, ||n(c)) is not orthogonalized but normalized, 
that is, {m{c)\\n{c)) / 6mn for m / n but (n(c)||n(c)) = 1. Since (c|c) = 1 and |n)) 
and ||n(c)) are the normalized states, we have the relation 

oo 

J2nn;cf = l. (5-7) 

n=0 

Since our intrinsic system which we are interested in is given only in terms of the 
6-boson, we need to evaluate an expectation value for any operator Ob consisting of 
only (6,6*) with respect to |c) : 

oo 

{Ob)c = {cmc) = J2rin;cf{n{c)\\db\\n{c)) . (5-8) 

n=0 

This situation is strongly similar to that in the case of thermo field dynamics for- 
malism. The relation ( ^-81) tells us that the factor r{n; c)^ denotes the weight of the 
mixture of the pure state ||?t-(c)) in the statistically mixed state under consideration 
and the quantum mechanical average gives us the statistical average for the physical 
quantity Ob automatically. In §9, the above will be discussed again. 

From the above argument, we learned that r{n; c)^ denotes the statistical weight. 
Then, it is possible to introduce the entropy in the present formalism. First, we note 
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the following relation : 

oo 

\c) =riha;c)Y,Mc))0\n)) , (5-9) 

n=0 

Ua = a*a . (5-10) 

Here, r{ha;c) is obtained by replacing n in r{n;c) with the a-boson number oper- 
ator. With use of the operator r{ha; c), we define the operator S in the form 

5 = -In r(n,;c)2 . (5-11) 

Then, the expectation value of S for the state |c) is calculated as follows : 

oo 

iS)c = {c\S\c) = - ^ r(n; cf lnr(n; c)^ . (5-12) 

ra=0 

The relation ( ^-12 ) tells that {S)c is the entropy and, then, S denotes the entropy 
operator. With the help of r(n; c) shown in Eq. (5-4), S is given as 

S = {ha + 1) In [/2 - fia In \V\^ + \Vf\6\^ - In iV(n,; \6f) . (5-13) 

The expectation value of S for |c) becomes 

{S)c = UHnU^ - \Vf\n\V\'^ + Ql^Fp^ {\V\'^lnU^ - \V\'^ln\V\^) 

|^f).i^ + .... ,5.14) 

Here, we used the relation 

InNiha, \6\^) = \5\^ ■ ha - \\5\Hha + hl) + --- . (5-15) 

We can show that, compared with the second term on the right-hand side of Eq. 
(5-14), the third one is expected to be small in a rather wide region of \V\'^. Then, 
hereafter, we treat the following form for the entropy : 



{S)c = U^lnU^ - \V\^ln\V\^ + (^^iW^l j {\V\^lnU^ - \V\^ln\V\^) . (5-16) 

With the use of the entropy S and the intrinsic Hamiltonian K, we define the 
free energy F in the form 

F={c\K-^\c) . (5-17) 

Here, (3 denotes the inverse of the temperature, i.e., /5 = (kBT)^^. The intrinsic 
Hamiltonian K in the present case is expressed as 

K = Ko + fieb*b + n'^gb*'^P . (5-18) 
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(i) harmonic oscillator 



5 = . 



(5-19) 



(ii) pairing model 



i^n = 



e = 2(e - GS) 



g = G. 



(5-20) 



(iii) Lipkin model 



^ h 3hG f ^(^ fi 

= -he + - cjo - [e + G[S + 

3G f ^f^n 

^0 , g = ^ [e + G[ S + 



The expectation value of K for |c) is calculated in the form 



(K)c = Ko + he(\V\'^ + -\W\'^U^ ] + Tv'g 



2\vr + ( -\w\ 



(5-21) 



(5-22) 



Here, we neglect the term with the order (2|VFp/?i)^. Then, combining {K)c with 
{S)c shown in Eq. (|5-1(^ ), we can express F as a function of \V\^. The variation of 
F with respect to \V\'^ leads us to the following relation : 



|y|2exp[/3f(|F|2)] = [/^exp 



Here, <5(|yp) is defined as 

£{\V\^) = he + Ah^g 



\V\' + 



\W\ 



1 + ||M^|2 C/2 



(5-23) 



(5-24) 



Under the approximation e" 
we have 



-X 



1-X for X = {2\W\'^/h)/{l + 2\W\'^/h) ■ (l/C/2 



1 



The relation ( ^-25 ) gives us 



n\v\ 



1+ 2|p^|2 ef(|y|2)/3_i ■ 



/ [/2 + |y|2 2|^|2 

In ' " 



|l/|2 + |y|2 2|^|2 



(5-25) 



(5-26) 



The quantities \W\'^ and \V\'^ are given in Eqs. ( 3-30| ) and (4-57), respectively, and 
Tz is obtained as a function of t. Therefore, with the aid of the relation ( [5 -261 ), f3, 
i.e., the temperature of the system can be determined as a function of t. 

In order to make further discussion simpler, let us treat the case of A = 0, where 
A is given in Eq. (4-27), i.e., the case of harmonic oscillator. As is shown in Eq. 
(|54g| ), Ko = 0, e = w and 5 = 0, and £{\V\'^) in Eq. (|]2|) becomes 

£{\V\^) = hoj . (5-27) 
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Then, the relation (5-25) is reduced to 



\V\ 



With the relations (5-22) and ( 5-2^ ), we have 



(e = hui) 



(5- 



W\^+{1 + -\W 
1 



1^1 



1 



(5-29) 



The quantity 2\W\'^ /h is a constant of motion and at the low temperature limit, the 
second term on the right-hand side of Eq. (5-29) vanishes. Therefore, the first term 
represents the energy at the low temperature limit and the second the energy coming 
from the thermal fluctuations in the boson distribution. 

Next, let us investigate the time-evolution of the system under a special initial 
condition (/ = /, J = 0,xo = 0). Then, |yp in Eq. (|^) can be expressed as 



= [sinh(7t)]2 . 



(5-30) 



Here, we used the relation (4-44) with 70 = 7. At the initial time t = 0, \V\'^ = 
and in this case r{n; c)^ becomes r{n = 0; c)^ = 1 and r{n 7^ 0; c)^ = 0. This 
means that the system is in a pure state with n = 0. At the time t — > 00, ^ 00 
and in this case, all r{n; c)^ becomes equal to one another. This means that the 
present system is in the equal weight mixing. On the basis of the above fact, let 
us investigate the time evolution of the system. The energy {K)c in Eq. (5-29) is 
expressed as 

{K)c = {21 -h)uo + uj- 2/|yp ( = E{t) ) . 
Since \V\^ = at i = 0, E{t) is expressed as 



E{Q) = {21 - n)uj . 
On the other hand, E{t) can be also expressed in the following form : 



E{t) = J2nn;cf{n{c)\\Kt\\n{c)) , 
n=0 

'dLn+i{x) 



(5-31) 
(5-32) 



{n{c)\\Kf,\\n{c)) = huju — huj\6\' 



dx 



-\5\ 



2 Ln{-m 



(5-33) 
(5-34) 



Then, we have the relations 

(Contribution of the 1st term in Eq.(5-34)) 



n=0 



(5-35) 
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(Contribution of the 2nd term in Eq.(5-34)) 



n=0 



' dLn+i{x) ' 
dx 



LU{2I - fl) . 



(5-36) 



It is interesting to compare the above results with the initial value of the energy 
given in Eq. ( p-32| ). The ensemble average of the second term shown in Eq. (5-36) 
coincides with the initial energy 

( F32D - The initial state is in the pure state with 



n = and the energy of the system concentrates on this state. At the time t (> 0), 
the states for n 7^ exist in |c) and the energy is distributed to each state. This 
fact tells us that the energy in the pure state with n = at t = dissipates to other 
pure states at t under the mixing weight r{n;c)'^. However, totally, it conserves. 
Next, let us investigate the contribution of the first term. It is seen from the result 
in Eq. (5-35) that this energy originates in the thermal effect. Therefore, the energy 
is supplied from the "external environment" and by this energy the temperature of 
the system increases. More detail informations are obtained in Refs.|l^), ^) and |2^. 

§6. The pseudo su{l, l)-algebra in the sM(2)-spin system 

As was discussed in §5, the su{l, l)-algebraic approach to the description of 
thermal effect gaves us various quite interesting features. However, this approach 
contains an undesirable feature. The Hamiltonian for the su{l, l)-spin system does 
not represent the energy of the entire system. Thus, it may not be justified to assume 
that the additional degree of freedom expressed in terms of the boson (a, a*) can be 
specified as a physical object characterizing the external environment. Rather, it 
may be natural to regard the additional degree of freedom as a theoretical tool 
for describing dynamical motion such as the damped oscillation in the frame of 
isolated system. In this sense, a fictitious degree of freedom is introduced in the 
su{l, l)-algebraic approach, in spite of giving us interesting results. Then, it may be 
interesting to present, without introducing the fictitious degree of freedom, a possible 
approach which gives us the same results as those obtained in the su{l, l)-algebraic 
approach^In this section, we will present the basic theoretical framework for this 
problem. c§) 

Let us prepare a boson space constructed in terms of the boson operators (c, c*) 
and {d,d*). In this space, we introduce the following operators : 



r+= Vhc* ■ VWTnFd ■ , _ ■ Vhd , 

\Jfi + M*d 

r_= Vhd* ■ _ • V2r + hc*c ■ Vhc , 
yh + M*d 

To=T + hc*c. (6-1) 
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Here, T is a real parameter which is assumed to be 



T = ht, t = ^ , 1 , ^ , 2 , ■■■ . (6-2) 



The above three operators are rewritten in the form 



1 o / o o 

= S+ V2T+ S + So 



n+s- So 

1 



T„= \/2T+ S + So- S- 



n+ s - So 

ro=T+S + So ■ (6-3) 

o 

The operators (^±,0) compose the sii(2)-algebra and are defined as 

S+= hc*d , S-= hd*c , So= |(c*c - d*d) , (6-4) 
S=^{c*c + d*d) . (6-5) 

o 00 

It is interesting to see that (t± o) can be expressed in terms of {S±,o^S)- Further, 
we can prove the following relations : 

0*00*0 

r_ =T+ , To =To , (6-6) 
[t+, t_] = -2h To +[(2T + fic*c){n + nc*c D) + {1- e)h ■ Titc] D , 

[to,t±] = ±nT± , (6-7) 

o2 1.00 00. 

^0 -2(^-'^+ + ^+'^-) 
= T{T-en) 

+]^[{2T + nc*c){n + nc*c)-{l-e)[2T + nc*c){n-nc*c)]D . (6-8) 

o 

Here, e and D are given as 

e = 1 , (6-9) 

o ^ I ^ o^o o*o 

D=l-d* ^.d, D =D , D =D . (6-10) 

l + d*d 



The relations (6-6) ~ (6-8) should be compared with the relations (3-2) ~ (|3-4| ). 

o 

If the terms related with the projection operator D are omitted, all the relations 
are reduced to those governed by the su(l, l)-algebra in the su(l, l)-space with 

o 

the magnitude T. The operator D plays a role of the projection operator, i.e., 
b •(l/\/n!) • ((i*)"|0))rf = Sno\0))d for d\0))d = 0. This means that the behavior 
of the set {t±^o) is quite similar to that of the su{l, l)-spin with the magnitude T. 



26 A. Kuriyama, J. da Providencia, Y. Tsue and M. Yamamura 

In this sense, it may be possible to regard the set (t±,o) as composing the pseudo 
su(l, l)-algebra. It is interesting to see that it is constructed in the su(2)-spin. 

Our next problem is to introduce a possible form of coherent state which leads 
us to the classical counterpart of the set (t-i-^o)- For the preparation, let us consider 
a state \a) which satisfies 

a\a) = a\a) . (6'11) 
Here, the operators a and a* are defined by 

a = f{a*a)~^a , a* = a*f{a*a)-'^ . (6-12) 

The operators a and a* denote boson annihilation and creation operators, respec- 
tively. Of course, f{a*a) and its inverse f{a*a)~^ should be defined. The eigenvalue 
a is, in general, complex. We can see that if f{a*a) = 1, the solution of Eq. ( |6T1| ) 
is the well-known boson coherent state. In association with a and a*, we introduce 
the operators A and A* in the form 

A = f{a*a)a , A* = a*f{a*a) . (6-13) 

The commutation relation between a and A* is given as 

[a,A*] = l. (6-14) 

Then, the state \a) is obtained in the following form : 

\a) = Mcexp{aA*)\0) , a|0) = , (a|a) = 1 . (6-15) 

oo 1 n— 1 

(M,)-2 = l + ^-n/(fc)'-(Hy . (6-16) 

n=l k=0 

For the complex variable z, we can derive the relation 

(alexia) = -{a*dza — adza*){a\f{a*a)\a) 

= ^{a*dza-adza*) . (6-17) 

Here, a is defined by 

a = a^{a\f{a*a)\a) . (6-18) 

The relation (6-17) supports that the parameters a and a* can be regarded as canon- 
ical variables in classical mechanics. The expectation value of a* a for \a) is given 
as 

(a|d*a|a) = {a\a* f {a* a)'^ a\a) 

= \a\'^{a\f{a*af\a) = a*a . (6-19) 

Under the above preparation, let us investigate a possible form of coherent state 
which gives us the classical counterpart of the set {T±fi). For this aim, we set up the 
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following wave packet : 



|c))=Meexp W)) , (6-20) 



\5)) = Jl-\6\^eKp 6A |0)) , (6-21) 



A =d*\ll + d*d . (6-22) 

Here, 7 and 5 denote the complex parameters. Clearly, the state \5)) is nothing but 
a simple example of \a) shown in Eq. (6-15) and satisfies 

5\5)) = 5\5)) , 5* = d* (^\jl + d*d^ , (6-23) 

T_ \6)) = . (6-24) 

Therefore, we can see that the state |c)) consists of a superposition of successive 
operation of t_|_ on the state \5)) satisfying the relation (6-24). With the use of the 
relation (6-23), together with the definition of shown in Eq. (6-1), the state |c)) 
can be rewritten in the form 



|c)) = (1 - |7p)*(l - . exp (7cV2t + c*c) exp y6(t^l + d*dj |0)) . (6-25) 

The relation (|648|) with /(c*c) = 2t + c*c and f{(td) = l + d*d tells us that 7 and 
6 can be expressed in terms of the canonical variables c and d : 



d 



(6-26) 



y/2t + c*c ' ^Jl + d*d 

The relation (6T9) gives us 

{{c\c*c\c)) = c*c , {{c\d*d\c)) = d*d . (6-27) 

O 

Then, the expectation value of S in Eq. (6-5) is calculated as 

{S)c = ((c| S \c)) = ^{c*c + d*d) = S . (6-28) 

Instead of the variables (c, c*) and (d, d*), we use the variables S) and {X, X* 
used in §2 : 

^ f ■'^\ * ^* f ■'4' 

c = X exp —i— I , c = X exp +i— 



d = -^V2S -nX*Xexp , d* = -^^28 - hX*X exp 



(6-29) 
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With the use of the relation (|]T|) for /(c*c) = V2t + c*c and f{d*d) = \ll + d*d 



together with the relation ( |6-26| ) and (6-29), the expectation values of (t±,o) for |c)) 
are given in the following form : 

' ^ ^ Vhx*V2T + nx*x ■ V2S - nx*x , 



Vh + 2S -hx*x 
(r_)c = V25 - hx*x ■ V2r + hx*x ■ Vhx- ^ 



y/h + 2S -hx*x ' 

(ro)c = hX*X + T . (6-30) 

We can prove that {T±fi)c obey the relations 

(r_): = (r+), , (ro): = (ro)^ , (6-31) 
[(r+)„ (r_),]p = i-i) {-2n(ro)c + [(2r + hc*c)ih + nc*cD) + (1 - e)?l • nc*c]D} , 
[{ro)c, {r±)c]p = {-i){±h{T±),) , (6-32) 

(ro)e'-^((T-)c(r+)e + (r+).(r_)e) 
= T(T - e/i) 

+i[(2T + hc*c){n + 7ic*c) - (1 - e)(2T + ?ic*c)(n - ;ic*c)]Z) . (6-33) 

Here, e and D are given as 

e = , (6-34) 
D = (l + d*d)-^ = l-d* ^—d. (6-35) 

It is interesting to see that the relations (6-31) ~ (6-35) correspond to the relations 
(6-6) ~ (6-10) completely. The quantity e appears from the non-commutability of 
various operators, i.e., it is a kind of quantum effect. From the above correspon- 
dence, we obtain a conclusion that the expectation values (t±^o)c are the classical 

counterparts of (t±^o)- 

With the help of the relation (6-3), together with the form (2-17), the Holstein- 
Primakoff boson representation for the set (t±^o) is obtained easily : 



1 



• Vnx*J2T + hx*x ■ J2S-nx*x , 



n + 2S- hx*x 

1 



f_ = \/2S-nx*x ■\/2T + nx*x ■Vnx^= 

\jh + 2S -nx*x 

fo = nX*X + T . (6-36) 

Here, S±fl and S in Eq. (6-3) are replaced with ^±,0 in Eq- (2-17) and 5, respectively. 
We can see that under the replacements {X^X*) {X,X*) and S S, the form 
(6-36) is reduced to the form (6-30). 
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§7. Description of a many-boson system interacting with an external 

harmonic oscillator 

The aim of this section was aheady mentioned in the beginning part of the pre- 
vious section. Corresponding the_FIamiltonian ( |5-18| ), we treat the following Hamil- 
tonian for a many-boson system E3): 



Kc= he"" ■ c*c + h^g^ ■ c*^c^ . 



(7-1) 



It may be self-evident why we treat the Hamiltonian (7-1) in the present case. The 
above boson system interacts with a harmonic oscillator, the Hamiltonian of which 
is given by 

Kd= fiuj^ • d*d ■ (7-2) 
As the interaction between the two systems, we are interested in the following form 



Vcd=Vcd (1)+ Vcd (2) , 

V^d (1) = fi^f ■ {c*d ■ d*c + (tc ■ c*d) , 

Vcd (2) = -?i7° • i 



h + hd*dj ■ c*d ■ V2T + hc*c 
V2T + hc*c ■ d*c ■ (^\/h + hd*d^ 



(7-3) 
(7-4) 



(7-5) 



Of course, e^, g^, to^, f^, and T denote real parameters. In particular, we regard 
T as the parameter obeying the condition T = ht, {t = 1/2, 1, 3/2, 2, • • •). Thus, the 
Hamiltonian of the entire system is expressed as 



Hsu(2)=Kc + Kd + Vcd 



(7-6) 



The meaning of the "sii(2)" in Hsu{2) is interpreted as follows : We are now describing 

o 

the system obeying the sn(2)-algebra and Hsu{2) can be expressed in terms of the 

o o , 

operators (£"±,0) a-nd S shown in Eqs. (6-4) and (6-5). The Hamiltonian ( [7- 61) can 
be rewritten in the form 



Hsu{2) = Usu{2) + Ksu{2) + Vsu{2) 



Usu{2) 



S 



-[(e°-a;0)-%°].T+(5°-2/0).r2 



(7-7) 



(7-8) 



ksu(2) = [(e° - - %° + 4/0 S -2(5° - 2f)T] To +(5° - 2/0). , (7-9) 
Vsu(2) = -l'' -iih-r^) ■ (7-10) 
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The operator 5* commutes with Hsu{2)^ ^iid U su{2) contains only S- Therefore, the 

o o 

dynamics of our present system is described in the framework of {Ksu(2) + V su{2)) 

o 

under a given eigenvalue 5 of S". 

In associating with the sti(2)-spin system introduced in the above, we will reca- 
pitulate a possible form of the su{l, l)-spin system discussed in §4. Our interest is 
in the following Hamiltonian : 

Hsu(l,l)=Kb-ka + Vba , (7-11) 

Kb = he- b*b + h'^g ■ b*H'^ , (7-12) 
Ka = he-a*a + h'^g-a*'^a'^ , (7-13) 
Vba = -h-/ ■ i{a*b* - ba) . (7-14) 

The Hamiltonian is of the same form as that given in Eq. ( |5-18| ) . The Hamiltonian 
(7-11) can be rewritten in the form 



Hsu(i,i) — Usu(i,i) + -^s«{i,i) + ) (7-15) 

2, 



Usu(i,i) = 2ie-hg){T-^) , (7-16) 



Ksu(i,i)=^9-{f-^^fo, (7-17) 

Vsu(iA) = -l-i{f+-f^) . (7-18) 

In a sense similar to the case of the SM(2)-spin system, the operator T commutes with 
Hsu(i,i): s-iid contains only T. Therefore, the dynamics of the present system 

is described in the framework of {Ksu{i,i) + ^su(i,i)) under a given eigenvalue T of 
T. As was already mentioned, the Hamiltonian Hsu(i,i) gives us various interesting 
aspects of thermal properties and dissipation for boson systems. In this sense, the 
form (7T5) is quite interesting. However, it does not represent the total energy of 

o ^ 

the system. On the other hand, Hsu{2) in Eq. ( [7 -61 ) represents the total energy of the 



system. Therefore, it is interesting to investigate the relation between both systems. 

At the present stage, we have two spin systems. One is the su{2)- and the other 
is the su{l, l)-spin system. For a given value of 5 (= hs; s = 0, 1/2, 1, 3/2, • • •), the 
STi(2)-spin system is described by the complete set {|n; 2s — n))}, which is expressed 
as 

In; 2s-n)) = ^_=L_(c*)"(d*)2^-"|0)) . (7-19) 
Here, n runs in the region 

n = 0, 1, 2, 2s . (7-20) 

For a given value of T (= ht, t = 1/2, 1, 3/2, 2, • • •), the su{l, l)-spin system can be 
treated by the complete set {|n; 2t — 1 + n)}, which is given by 

|n; 2t - 1 + n) = , , ^ (a* f ~^+"|0) . (7-21) 

^/n\[2t — 1 + nj! 
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Here, n runs in the region 

n = 0, 1, 2, 2s, 2s + l,--- . (7-22) 



We refer to the spaces spanned by the complete set ( 7-19|) and ( 7-21| ) as the su{2)- 



space with s and the su{l, l)-space with t, respectively. 

Let us transcribe the su(2)-spin system of the sii(2)-space with s into that in 
the su{l, l)-space with t. For this aim, first, we set up the following correspondence 

|n;2s-n)) ~ |n;2t- 1 + , (7-23) 
n = 0, 1,2, •••,2s . (7^24) 

Following the basic idea of the boson mapping, we call the space spanned by the set 
{|n; 2t — 1 + n);n = 0, 1, 2, • • • , 2s} as the physical space and define the operator Ugt 

2s 

U,t = J2\n;2t-l + n){{n;2s-n\. {7-25) 
n=0 

The operator Ust satisfies 

2s 



ulAt = E 2s - n))((n; 2s - n| = 1 , (7^26) 

n=0 
2s 

UstUlt = ^ |n; 2t - 1 + n){n; 2t - 1 + n\ = Pst ■ {7-27) 



n=Q 

Here, Pst plays the role of the projection operator to the physical space {Pg^ = 
Pst, Pst = Pst)- The operators Ust and JJI^ give us 

iJst\n]2s-n)) = \n;2t-l + n) , (7^28) 

U\\n-2t-l + n) -l^'^'''^^''^^^ ' (n = 0,l,2,^^^2s) 

U,t\n,lt i + n)-<^^^ {n = 2s + l,2s + 2,---) ^' 

o o 

The operator Ust permits us to transform t±^q, S and T as 

Usf h,o -Ult = Pst ■ T±,o • Pst , 

Usf S -ult = UsfS- Ul = PsfS- Pst , 

UsfT- ill = PsfT- Pst = Psff- Pst . (7^30) 

Here, it should be noted that S and T are c-numbers. With the aid of Ust aud the 

o 

relation (7^30), we are able to transcribe Hsu{2) iu the sm(1, l)-space with t : 

^ O - J. 

Usf Hsu{2) -Ust = Pst ■ Hsu{2) ■ Pst , (7-31) 
Hsu{2) = Usu{2) + Ksu{2) + ysu{2) , {7-32) 
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■S 



_[(eO _ _ ^^0] . f + (5° - 2/0) . f 2 , (7-33) 
Ksui2) = [(e° - u') - hg' + AfS - 2{g^ - 2f)f]fo + (5° - 2/°) • f ^ , (7-34) 
V;«(2) = -7'-^(r+-r-). (7-35) 

We can regard -^3^(2) the Hamiltonian in the l)-spacc which has a connection 

o 

with the original Hsu{2) i^i sti(2)-space. In the same meaning as that mentioned 
for the Hamiltonian ^^3^(1,1)) we are interested in the part {Ksu{2) + ^su(2))- Then, 
we require the following condition for the parameters {e^,uj^,g^, f^,j^,g,j) : 

e»-.» + 2,»(5-i)=4,(T-f), 
g° = 2f\ 7" =7. (7-36) 
Under the condition (7-36), we obtain 

Pst ■ {Ksu(2) + Kw(2)) • Pst = Pst ■ {Ksu{l,l) + Kw(l,l)) ■ Pst ■ (7-37) 

This means that, under the condition (7-36), Hsu{2) is expected to display various 
results similar to those given by Hsu{i,i)- ^ possible reexpression of the condition 
(7-36) is given as follows : 

^0 _ 50 frr. ^\ n frr ^ 



7° =7, 9 = f ■(■S-f) ■ (»0 : parameter) (7-38) 

o o ^ ^ 

Under the condition (7-38), {Ksu{2)}Vsu{2)) a-^d ^su(i,i)) can be expressed 

in the form 



Vsu{2)=—fiir+-r^), (7-39) 

Ksuii,i)=2go{f-^) (5-^) -To, 

Vsu{i,i) = -l-iiT+-f_) . (7-40) 

Thus, in the s'u(2)-space with S and the 5^(1, l)-space with T, we have the following 
correspondence 

o ,^ o ~ 

Ksui2)^ Ksu{l,l) , Vsu{2)^ • (7-41) 

The above correspondence results from 

r±,o~ r±,o , T ~ f , S^S . (7-42) 
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Of course, Usu(2) and C/^„(i,i) 



noted that there exists a difference between Hsu{2) and Hsu{i,i) 



do not correspond in forms. However, it should be 

o 

in addition to U su{2) 

and tJsu{i,i)- The Holstein-Primakoff representations shown in Eqs. (3-16) and (6-36) 
teach us the difference. In the case of (7±^o) in Eq. (3-16), the integer n in the boson 
state (l/-v/nl) • (X*)"|0) is permissible up to the infinity. But, in the case of (t±^o) in 



Eq. (6-36), n is permissible up to 2s. This can be seen in the factor y2S — hX*X. 
It is quite natural, because the sn(2)-algebra is compact, and the su{l, l)-algebra 
noncompact. If s is sufficiently large compared with n, the two algebras give almost 
the same results. 

Next, we investigate the correspondence between the two systems in the classi- 

o 

cal counterpart. The expectation values of Hsu{2) ^^'^ Hsu{i,i) with respect to the 



coherent states (|6-25| ) and (|3-23|) , respectively, give us the following results : 



su(2) 



u. 



su{2) 



su{2) -r Vsu{2) ; 



U. 



su(2) 



2f. 



T -- 
2 



■ S 



5^1 



K. 



su{2) 



%° + 4/05-2(50-2/0)r](ro)c 



5^1 + 



n 

2T 



2f 



Vsui2) = -l°-ii{r+)c-{T-)c) 



H 



su{l,l) 



u. 



u. 



su{l,l) 



su{l,l) + -^su{l,l) + ) 

2{e-2hg) (t-^^ 



45 



h 



h 



77 1 + :^ (^C 



2T 



'0)c 



(7-43) 



(7-44) 



(7-45) 

(7-46) 
(7-47) 

(7-48) 

(7-49) 
(7-50) 



Here, (t-i-^o)c and (7±^o)c are given in the relations (6-30) and (3-32), respectively. In 
the same meaning as that in the quantum version, Usu{2) and Usu{i,i) do not relate 
with the dynamics. For the Hamiltonian (i^^5u(2) + ^su(2)) and (^^^^(i,!) + 
we set up the following condition : 



UJ 



90 



/0 = 5o /^_-| /^^^ 



T 



h 



2T 



90 
2 



S . 



(7-51) 



The relation (7-51) corresponds to Eq. (7-38). Then, {Kg^2)^ ^su(2)) and (i^sn{i,i); ^sn{i 
can be expressed in the form 



su{2) 



h 



h 



25o(r--) + — )5(ro) 



(7-52) 
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Vsu{2) = -7 • «((t+)c - (t_)c) , (7-53) 

Ksu(i,i) = 250 {t - I) (l + ^) S{n)c , (7-54) 

Kn(i,i) = -7-^((r+)c-(T_),). (7-55) 

It should be noted that (S*, T) are the parameters and {S^ T) denote the dynamical 
variables, which are constants of motion in our present system. Under the condition 
{S = S,T = T) and the correspondence (t-i-^o)c ~ (7±,o)c, we have 

Ksu{2) ~ -f^sM(l,l) ) ^sn(2) ~ Ku{l,l) ■ (7-56) 

However, in the case of the su(2)-spin system, we have the restriction < 2S 

and in the case of the su{l, l)-spin system, these does not exist such a restriction. 
Then, if S is sufficiently large compared with hX*X, the results of the two systems 
are almost the same as each other. 

Let us investigate the difference between (r±)c in Eq. (6-30) and {T±)c in Eq. 
(3-32). Under T = T, both are related with each other through the relation 



e{hx*x) 



(r±)c = {T±)c^l - e{hX*X) , (7-57) 

= 2S + n-nx*x • ^'-''^ 

The behavior of the function 0{hX*X) is summarized as follows : 

0{h+^) , {nx*x < 25) 
o(7i+5) , {hx*x ^2S -0{h+^)) 

0{h^) , {hX*X ^2S ±0{h+^)) (7-59) 

0{h^^2) , {hx*x ^2S + h-0{h+l)) 
[o{h-^). {hx*x ^2S + h-0{h+^)) 

Here, 0{h^) denotes the term of order h^. Therefore, under the condition 7i — > 0, 
6{hX*X) can be approximated in the form 

r , {hX*X < 25) 

e{nx*x) = < 1 , {nx*x = 25) (7-60) 
i 00 . lnx*x > 25) 

The behavior of 9{TiX*X) given in the above tells us that 6{hX*X) is equal to zero 
in the region hX*X < 25 and at the point fiX*X = 25, 6{hX*X) behaves like an 
infinite wall. The existence of the wall suggests us that the quantity hX*X shows a 
certain periodic behavior with the period ~ In 5. Then, if 5 is sufficiently large, we 
can expect that the solution for the pseudo su{l, l)-algebra is well approximated as 
that for the su{l, l)-algebra. The details are discussed in Ref. |2F 



§8. The su{2, l)-algebra and its modified form in terms of three kinds of 

boson operators 



In the previous section, we presented a possible description of a boson system 
interacting with an external harmonic oscillator. Of course, two kinds of boson 
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operators (c, c*) and {d,d*) are introduced and a coherent state for a possible time- 
evolution of the system is expressed in terms of these two kinds of boson operators, 
which describe the physical objects. Therefore, these bosons are not related to the 
phase space doubling such as mentioned in §5. In this sense, with the help of the 
coherent state adopted in §7, it may be impossible to treat statistically mixed state. 
In this section, we will give a theoretical framework which enables us to describe the 
mixed state of a boson system interacting with the external harmonic oscillator. IHs^ 
In addition to the boson operators (c, c*) and {d,d*), we introduce the boson 
(a, a*). With the use of these bosons, the following bi-linear form can be defined : 

T+= hac* , T-= hcd , fo= |(a*a + c*c) + | , (8-1) 

R+= hd*d* , R^= Ma , Ro= ^(a*a + d*d) + | , (8-2) 

(T+= hc*d , a-_= hd*c , (Jo= ■^(c*c — d*d) , (8-3) 

o o o 

The set (r±,o) and (ii±,o) obey the sii(l, l)-algebra, respectively, and (o"±,o) the 
sit(2)-algebra. The operator ctq is related with 

ao=To - Ro ■ (8-4) 

O O o 

We can prove that the set {T±fl, R±fl, cr±) composes the sm(2, l)-algebra. In the 

O O o 

above set, there exists an operator which commutes with T±,0) R±,o and a±^Q : 
X=^{c*c + d*d-d*a) + h , (8-5) 
[K, T±,o] = [K, R±,o] = [K, a±,o] = . (8-6) 



As was shown in Ref .|28|) , there exist many cases for constructing the orthogonal 
set of the su{2, l)-algebra. In this paper, we show one example. For this aim, we 

o 

introduce the operator T in the form 

T=^ic*c-d*a) + ^. (8-7) 

o o o 

The operator T commutes with K and To and, then, we set up the eigenvalue 
equation 

K \K,T,To)) = K\K,T,To)) , K = hk , 
T\K,T,To)) = T\K,T,n)) , T = fit , 
To \K,T,To)) = To\K,T,To)) , Tq = ^to • (8-8) 
A possible solution of the eigenvalue equation (8-8) is given as follows : 

\K,T,To)) ^ 



y/{to + t-iy.i2k-2t-iy.{tQ-ty. 
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x(c*)*«+*-^(d*)2'=-2*-^(a*)*o-*|0)) 
1 1 (2/c-2t- 1)! 



(2fc-2)!V (io+t-l)!(to-t)! 

(o \ io-t. / o \ 2i-l 

l^j {d*r-'m. (8-9) 

k = 1, — , 2, • • • , 

_ 1 3 k-l 
*~ 2' 2''"' ' 

to = t, t + t + 2,---. (8-10) 

The classical counterpart can be derived by the following coherent state : 

' o \ / o ^ 

|c)) = M.exp I ^ • ^ exp L)^ . ^ I expM*)|0)) . (8-11) 

Here, V , w and v are complex parameters and U = ^/T+'\V\^. Instead of the above 
parameters, {6,K), {X,X*) and (Y,Y*) are convenient for our discussion. They 
are required to satisfy the relation 

{{c\ihdQ\c)) = K - h , {{c\indj^\c)) = , (8-12) 

{{c\dx\c)) = +^ , {{c\dx'\c)) = -f , (8-13) 

{{c\dY\c)) = +^ , {{c\dY'\c)) = -y . (8-14) 

Further, K obeys the condition 

((c| k |c)) = K . (8-15) 
The relations (8-12) ~ ( ^-15 ) give us the form 



w = exp [ —I— ) A 



/ e\ 2{K-h) 

v = e.p (,-^2jV^ ' 

V = . (8-16) 

Vl + X*X 



Then, the expectation values of T±,o and R±fl are given by 

(r+)c = VhY* ■ Vhx*x + n + nY*Y , 
(r_)c = \/hx*x + h + hY*Y ■ Vky , 
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{To)c = hY*Y + ^ihX*X + h) 



= VhY* ■ p{K -h)- nx*x . ^^^.^ ^ ^ • 



{R.)c = Vnx* ■ ^^^.^ ^ ^ • \/2(K -h)- hx*x ■ Vhv , 

{Ro), = K-^X*X + hY*Y-^. (8-17) 

In Ref -p^) , it was discussed in detail that the form (8-17) can be regarded as the 
classical counterpart of the su{2, l)-algebra in terms of three kinds of boson opera- 
tors. 

In order to describe a many-boson system interacting with an external harmonic 
oscillator, we modified the su{l, l)-algebra, which was called the pseudo su{l, 1)- 
algebra. Then, in the present case, we must modify the su{2, l)-algebra. For this 
aim, first, we replace <t±,o in Eq. (8-3) with T±fi defined in the relation (6-1). Then, 
by extending the coherent states shown in Eqs. (6-20) and (|8-11 ), respectively, we 
adopt the following form : 

|c)) = M.exp I ^ • exp (^-^ ■ exp (^r? |0)) . (8-18) 

o * 

Here, V, ^ and r/ are complex parameters and U = ^yl -\- \ V\'^. The operator A 
is defined in Eq.(6-22). In the same way as in the previous case, we introduce 



{0,K), {X,X*) and {Y,Y*) which obey the conditions (8-12) ~ (|8T5D . Then, the 
parameters ^, r] and V are expressed in the form 

4 = exp —I— X - 



2 J V2t + X*X 



e 



r? = exp -i- J2{K - h) - hX*X 



2/ V ^ ' ^2{K -h) + n-hx*x 

V = Y , ^ . (8-19) 

On the basis of the form presented in the above, we will treat statistically mixed 
state of a boson system interacting with an external harmonic oscillator. The state 
|c)) in Eq. (pT8|) can be factorized in the form 



|c)) = l^-ca)) ® Wd)) , (8-20) 

^d^A . exp (i< 



|^ca)) = t/"'(l-|el')*exp(^G*c*') •expf|cV2t + c*c) |0))ea , (8-21) 



^d)) = (1 - \r]VY^ exp r/(i*Vl + d*d |0))d . (8-22) 



The above factorization will be used in the next section. 
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§9. Description of statistically mixed state of a boson system 
interacting with an external harmonic oscillator 

Before entering the concrete treatment for the present problem, we will contact 
with the basic idea of our treatment in rather general form. Let us investigate a 
many-body system consisting of two parts. We refer to the first and the second part 
as the c-part and the d-part. The Hamiltonian consists of 

H=kc + kd + Vcd ■ (9-1) 

o o 

Here, Kc and Kd denote the Hamiltonian of the c- and the d-part, respectively, and 

o 

the two parts interact with each other through the term Vcd- Therefore, we must 
introduce two spaces, which we call the c- and the d-space, and the present system 

o o 

is treated in the product of the two spaces. The operators Kc and Kd belong to the 

o 

c- and the d-space, respectively, and Vcd connects the two spaces. 

For the Hamiltonian (|9-1|), we set up the time-dependent Schrodinger equation 



indt\c)) =H \c)) , ((c|c)) = 1 . (9-2) 
As a possible solution of Eq. (|9-2|), let us assume the form 



|c)) = |<P,)) l^',)) , 

{^cp,)) = ii^dWd)) = 1 . (9-3) 

Here, \'Pc)) and \^d)) are given in the c- and the d-space, respectively. For the state 
(9-3), we impose two conditions. The first is the following : For any state \k)) in the 
d-space which is orthogonal to \^d)), \^d)) is governed by 

m kd + Vcd -ifidtWd)) = , (9-4) 
{{kWd)) = . (9-5) 

The second condition is the following : For the state \^d)) satisfying the condition 
(9-4), 1^,)) obeys 

indt\<Pc)) =Hc |<^c)) , (9-6) 
Hc=kc kd + Vcd -iMtl^d)) ■ (9-7) 

Here, it should be noted that the left-hand side of Eq. (9-4) and the second term 
on the right-hand side of Eq. (9-7) are calculated only for operators in the d-space 
and, then, naturally, they are operators in the c-space. Under the conditions (9-4) 
and (9-6), we can prove that the state |c)) given in Eq. (9-3) is an exact solution of 
the Schrodinger equation ( |9-2| ). As is clear from the forms (9-4) ~ (9-7), the above- 
mentioned two conditions are too strong to get the exact solution. However, it 
may be permissible as an approximation. In general, the time-dependent variational 
procedure gives us a possible approximate solution of the Schrodinger equation. 
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Therefore, under the procedure with the trial function |c)) in the form (9-3), we 
obtain an approximate solution of Eq. (|9-2D . 

Our chief concern is related to a possible description of statistically mixed state 
realized in the c-part. Then, with an idea similar to that used in thermo field 
dynamics formalism, we perform the phase space doubling for the c-part. This 
procedure is done through the introduction of a new space which is conjugate to the 
c-space. We call this space the a-space. Let a possible solution of the Schrodinger 
equation in a form analogous to |c)) given in Eq. (9-3) be obtained : 

\c)) = \^^))0\^d)) , {{c\c)) = l. (9-8) 

The part jlP'd)) obeys the conditions (9-4) and (9-5) and |^ca)) is given by solving 
Eq. (9-6) in the ca-space : 

ihdtl^ca)) =Hc l^ca)) , ii^cal^ca)) = 1 • (9-9) 

o 

Here, He is given in Eq.(9-7). With the use of the state |c)) obtained in the above 
process, we are able to describe the mixed state in a manner similar to the stan- 
dard thermo field dynamics formalism. In our present system, let us introduce the 

o 

operator Pc in the form 

^^=^(n|^,,))((^,,|n) . (9-10) 

n 

o 

Here, {|n)} denotes an orthonormal set in the o-space. Then, is an operator given 

o 

in the c-space. We can prove that Pc satisfies the conditions which the density matrix 
obeys : 

Tr h= 1 , ifidt Pc= [Hc,k] ■ (9-11) 

o o 

Here, He is given in Eq.(9-7). The ensemble average of the operator Oe in the c-space 
is expressed as 

Tr (dePe) = {{^ea\ Oe l^ca)) = ((c| Oe |c)) . (9-12) 
o 

Further, the entropy Sc in the c-space can be defined as 

Sc= - J2{{^ea\n){n\^ea)) ' ln( (^,, |n) (n|^ea) ) . (9-13) 

n 

We can discriminate the pure state in terms of the index n. 

On the basis of the basic idea mentioned in the above, let us investigate the 
Hamiltonian ( [7-6D in the framework of the modification of the su{2, l)-algebra given 

o o o 

in §8. The expectation values of Ke, Kd and Ved (1) shown in Eqs. ( |7-1| ) ~ (7-3), 
respectively, are calculated as follows : 

{Ke)e = e° • h{X*X + Y*Y) 

+g° • |2[n(X*X + Y*Y)f - - ^) {hX*Xf[l + hY*Y{h + hx*xf]j , 
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(9-14) 

{Kd)c = • [2{K -%)- nX*X] , (9-15) 
(Krf(l))c = 2/° • n{X*X + Y*Y)[2{K -h)- hx*x] + hf ■ [2{K -h) + nY*Y] . 

(9-16) 

The above three are exact. But, in the case of {ycd{2))c-, the approximation form is 
given as 



/ 2{K-n)-hx*x 
2{K -h) + h-hx*x 



(9-17) 



{Vcd{2))c = -7° • iVn{X* - X)V2T + hX*X^ 
Here, r is defined by 

T = T + nY*Y, T = nt. (9-18) 

Concerning the derivation of the relation ( |9-17| ), the following approximate form 
plays a central role : 



{{c\VhcW2T + hc*c\c)) = exp fi-^ ■ VhXW2T + hX*X 



(9-19) 



The detail is given in Ref.29). Further, we introduce the quantities T±fi as follows : 



r+ = VhX*V2T + hX*X^ 



I 2{K -fi)- hX*X 

2{K -h) + n-hx*x 



rrr i , , , / 2(K-K)-'hX*X 



To = r + hX*X , 

S = K-h+-Y*Y . 
2 



2{K - fi) + n - nx* X 



(9-20) 
(9-21) 



o ^ ^ 

The quantity S is the expectation value of S (= fi{c*c + d*d)/2) for the state |c)). 
Then, we have 



Hsu(2,l) — Usu{2,l) + -f^su(2,l) + ^sm(2,1) 



u 



su{2,l) 



= 2 



S 



-hg'^].T + 



Ksui2,i) = [(e° - - %° + 4/O5 - 2(5° - 2f)T]To 



+ 



1 + ^1-2/" 



'0 > 



(9-22) 

(9-23) 

(9-24) 
(9-25) 



We can see that the Hamiltonian (9-22) with the forms (9-23) ~ (9-25) is of the same 
form as that shown in Eq. (7-43) with the forms (7-44) ~ (7-46) if r and {K — K) in 
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r-t.o appearing in Eq. (9-20) are replaced with T and S in (t-i-^o)c appearing in Eq. 
(6-30), respectively. However, as was discussed in Ref J2^) , the meaning of Hgu{2,i) 
is different from that of /?s«(2)- Since ^^^su(2,i) is of the same form as -ff<ju(2)) '^^ can 
treat -ffsu(2,i) under the same method as that adopted in §§4 and 7. Then, we will 
not contact with this method in detail. 

In order to make the problem simpler, we consider the case of e^ = u:^ and 
gO _ jO _ which corresponds to the harmonic oscillator. Further, {K — h) is 
sufficiently large, i.e., 2{K — h) > hX*X. In this case, hX*X shows a certain long 



periodic behavior and, with the help of Eqs. (4-44), (4-47) and (4-56), we have 



hX*X = ^I' + cosh(27t + X 



(9-26) 



I.e., 



{{c\hc*c\c)) = + cosh(27i + x°) - T . (9-27) 
Since 6 = 2ujt + 6^, to order 7^, the expectation value {{c\hc*c\c)) can be expressed 



as 



{{c\hc*c\c)) = {{c\hc*c\c))t=o 



—ft 



{{c\Vhc* ■ V2T + nc*c ■ exp 



+V2T + hc*c ■ Vhc ■ exp ( i— j |c)) 



. (9-28) 



t=o 



Here, the symbol t denotes time. 

Now let us investigate our system in the framework mentioned in the beginning 
part of this section. With the use of \^ca)) and jlZ'd)) in Eqs. (8-21) and (8-22) and 

o 

the form r/ in Eq. (8-19), He defined in the relation (9-7) can be expressed as 



Hc = Ho + Hi+H2 , 

o 

Ho = hujc*c , 

7 • i[Vhc* ■ V2T + hc*c ■ 7] - V2T + hc*c ■ Vhc ■ r/*] 
+ i(jjrj)rf — [f]* — iujri*)r]] . 



Hi 

O 

H2 



ih 



1 



(9-29) 
(9-30) 
(9-31) 
(9-32) 



2 1- 

Substituting 6 = 2ojt + 6*° into the form (8-19), r] can be expressed as 

^ = (-4) " ^ - ^) - ■ 

p o ^ 

Then, we can prove that, at order h , the part H2 vanishes and rj = exp{—i9 /2) 

o 

X exp(— iwt). Thus, at this order. He can be approximated as 



(9-33) 



He = Ho + Hi{t) 



(9-34) 
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o o 

Here, Ho is given in Eq. (9-30) and Hi{t) is written in the form 



Hi{t) 



-7 • I 



Vhc* ■ \/2T + c*c ■ exp(-i6'°/2) • exp(-ia;t) 
-V2T + hc*c ■ Vhc ■ exp(+i6'°/2) • exp{+iujt) 



. (9-35) 



The Hamiltonian He indicates that the present system is a harmonic oscillator (or 
many-boson system) in a classical external field with the frequency uj. The inter- 

o 

action with this is Hi{t). We will apply linear response theory of non-equilibrium 
statistical mechanics. For the Hamiltonian (9-29), this method leads us to the fol- 
lowing form in the case of hc*c : 

TT{hc*c P {t)) = Tiihc*c °P (0)) 



-7t • Tr Vhc* ■ V'2T + hc*c ■ exp 



-V2T + hc*c ■ Vhc ■ exp +i 



Here, P {t) obeys the von Neumann equation 

ihdtPit) = [Hc°P{t)] 




(9-36) 



(9-37) 



o . . 

A possible solution of Eq. (9-37) is given in the form p (t) shown in Eq. ( 9-lOD and 
the relation (|9-12| ) tells us that, for any operator Oc, we have 

Tr(0, °p (0)) = {{c\d,\c))t=o ■ (9-38) 

Then, the results (9-28) and (9-36) coincide with each other under the correspondence 

(Q). In our present method, we can calculate Tr(7ic*c P (t)) for any order of 7 
in the specific order. Also, it may be interesting to see that our approach is closely 
related to linear response theory of non-equilibrium statistical mechanics. clP 

In our present system, hY*Y does not depend on t, while hX*X depends on t. 
Therefore, in the present system, the minimum value of hX*X for t can be regarded 
as the equilibrium value, that is, the equilibrium state appears. This state realizes 
at 27t + x° = in the relation (|9-26| ). Then, we have 



(hX*X' 



cq 



(T + hY*Y)^ + J^-{T + hY*Y) . 



(9-39) 



In this case, the entropy Sc is derived from the relation ( |9-13| ) in the following 
approximate form : 



Sc = {l + Y*Y) ln(l + 1^*1^) - Y*Y ln{Y*Y) 



(9-40) 
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At the equilibrium point, we define the free energy Fc in the form 

F, = E,-^ , (9-41) 
Ec = {{c\hLJc*c\c)) 



= uji^^JiT + nY*Y)^ + J'^ -Tj , (9-42) 

P = (fcBTeq)-l . (9-43) 

By minimizing Fc with respect to Y*Y, we have the following relation : 

I = T + n y ^ . (9-44) 



exp [hu^P . j - 1 

From the above relation, we can calculate the energy Ec as a function of huP. Two 
extreme cases are rather easily derived : 

(1) The low temperature limit : 

Ec = , „ + nujc 



T + Vt^TJ^ exp{huJcP) - 1 

h huJ^ / 1 Y 

"2(Vr2 + J^)3 ■ Uxp(?iu;,/3)-lj ' ^' 

ujc = uj ^ J ^ . (9-46) 



(2) The high temperature limit : 

Ec = kBT,q + uj^\It^ + J^ ■ {kBT,^)-^ + ■■■ . (9-47) 

We see that in the relation (9-45), the first term denotes the energy of T^q = 0, 
and the second represents the energy coming from the thermal fluctuations and the 
distribution function of free bose particle. The frequency w is changed effectively 
to LOc, which may be the result of interaction with the external oscillator. On the 
other hand, the high temperature limit reveals that Ec ~ k^T^q and we can see that 
present system obeys the well-known law of equi-partition of energy. The details can 
be found in RefJH). 



§10. Time-dependent variational approach with squeezed state to 

quantum mechanical systems 

In this and succeeding sections, we introduce another trial state in the time- 
dependent variational approach for quantum many-particle systems with the aim of 
investigating the classical motion including quantum fluctuations in the quantized 
systems. The new state, that is called the (one-mode) squeezed state, is introduced in 
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this section, which gives a possible classical motion in quantum mechanical systems 
including quantum fluctuations beyond the order of TiM3 The similar variational 
approach is found in Ref.^), while our treatment is developed widely in the system- 
atic way. We here present instructively our formalism based on the time-dependent 
variational principle with the squeezed state in quantum mechanical systems. 
The squeezed state is defined as 



\^{a,(5)) = (1 -/?*/?) texp( ^6*2 



• exp ( — a* a] exp(aa* 



(10-1) 



Here, a* is a boson creation operator and |0) is a vacuum state for the boson anni- 
hilation operator d : d|0) = 0. The operators h* and h are defined as 



a 



a — a . 



(10-2) 



The operator b is nothing but the annihilation operator for the usual coherent state : 
6exp (^—^a*a^ exp(ad*)|0) = 0. Introducing the operators which correspond to the 
coordinate and momentum operators, we have another expression of the squeezed 
state ( p4| ) : 



(a,/5)) = e'^(2G)-i/Sxp ( -{pQ - qP) ) exp 



Q 



'-X(f+a) 



1 

2^ 



2G 



P 



+ i277 ) Q2j |o) , 
(10-3) 

(10-4) 



G 



-{a* + a) 



+ \y\ +y 



-i2ip 



G 



P 



n 



a 



iy*-y)\ - + \y\^G~\ (10-5) 



- + |yr + y 



where y is related to /? as y = /3/-\/2(l — The reason why we have introduced 

the new variables y and y* instead of /3 and /?* is shown later. 

The time-evolution of this quantum state is governed by the time-dependent 
variational principle : 



to 



dt{i;{a,P)\indt- H\i^ia,l3)) = . 



(10-6) 



In the present method, one can deal with any Hamiltonian which is a function Q 
and P. However, we restrict ourselves in this section on the Hamiltonian of the form 



H=\p-^ + V{Q) . 



(10-7) 
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Further, the following canonicity conditions 



1 



(V(a,/3)|9,|^(a,/3)) = -x* 



{^l^ia,(3)\^y\^P{a,(3)) = -y* , 
give us the sets of canonical variables 



X = a 



a 



/3 



{iP{a,P)\d,.\i;{a,(3)) 
(V(a,/3)|5,*|V'(«,/5)> 

f3* 



1 

-r 
1 

-2^ 



V2(l -/?*/?) ' 



V2(l -/?*/?) 



(10-^ 



(10-9) 



This is the reason why we have already expressed G and U in Eq.(10-5) in terms of 
y and y* instead of /3 and /?*. Thus, the variational principle ( |10-6| ) leads us to the 
following equations of motion : 



q = dpH 



P ■ 



-dqH 



ihy = dy*Ti. , ihy* 



-dyU , 



(10-10) 



where 7i is the expectation value of the Hamiltonian with respect to the squeezed 
state. Also, the variables G and 77 obey the similar canonical equations of motion 
instead of the variables y and y* : hG = djj'H and hU = —dc'H. Here, we summarize 
the expectation values for the various operators : 



{^{a,(5)\Q\tlj{a,f3)) = q, 
(V(a,/3)|QV(«,/9)) = 9' + ^e 



(V'(a,/3)|P|V(a,/3)) = P, 
(V^(a,/3)|p2|^(a,/3)) = + n?? , (10-11) 



(V'(a,/3)|i/|^/'(a,/3)) = ttP^ + tt^?? + exp ^ -/i^ 



d_ 
dq 



V{q) 



(10-12) 



where the real quantities ^ (> 0) and ij {> 0) are defined as 



- + y*y + y 



G 



- + y*y-y 



^10-13) 



The squares of the standard deviations for Q and P are expressed as (^/^(a, /3)|(Q — 
g)^|V'(a, /3)) = and {'4}{a^ f3)\{P — p)^|V'(a, /?)) = hrj. Thus, the uncertainty re- 
lation is expressed in terms of ^ and rj as {ijj{a, I3)\{Q — q)'^\ijj{a, fi)) {ilj{a, [5)\{P — 
p)'^\il){oi, (3)) = h'^^rj = h"^ {1 / A + G"^ n^) . It can be seen from this uncertainty relation 
that, if one direction of the uncertainty is relaxed, the other can be squeezed. Also, 
one can see from (10-12) that the quantum effects beyond the order of h are included 
in this formalism. This novel feature is originated from the degree of freedom of the 
squeezing. 

We have to give the initial conditions for solving the Hamilton equations of 
motion derived in Eq. ( 10-10| ). Among four initial values for qo, pQ, yo and ?/g which 
lead to a unique set of solutions of Eq.( |10-10 ), those for qo and po can be selected 
arbitrary as is usually done in the ordinary TDHF theory. Thus, our task is to 
determine the initial conditions for y and y* . Since we are interested in classical 
motions in quantal systems, it may be natural to require that, at least, the initial 
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state has the least quantal effects. Thus, the initial state satisfies the minimal 
uncertainty relation. Also, this means that the absolute value of the energy due to 
quantum effect at the very beginning is as small as possible. These two conditions 
give us the initial values of yo and : 



1 



+ VoVo - Vo 



1 
4 



\'H{qo,Po,yo,yo) - Wc(go,Po)| : minimal . 

(10-14) 

Here, each variable with suffix denotes the initial value of each dynamical variable, 
and Ti-dq^p) denotes the classical value of Ti defined as 

1 



(10-15) 



The first condition in Eq.( |10-l"^ ) gives us yQ = yo (= p), which is identical with 
TTq = 0. The second condition in Eq. (10-14) together with y^ = yo is reduced to 



1 



1 

8Go ^ h 



exp { ^hGo { 1^ 



minimal. 



(10-16) 



Here, Go has already been defined in the relation (10-5). 

In the order of and h^, the condition ( |10-16| ) is further reduced to 



1 

Go 



+ 4vGo + hwGl 



minimal 



a2 

7rT^(9o) 



o4 

w = j-^Viqo). (10-17) 



Approximate solutions of Eq. (|10-17) and the corresponding energies are given as : 



i) ?; > 0; 




Go 


1 hw 
2^ 16^2 ' 


E = 


Ecfi + 


^ w 

2 62 V 


(10-18) 


ii) t; < 0; 




Go 


1 hw 

~ 2^ ' 32^;2 ' 


E = 






(10-19) 


iii) V = 


,w > 0; 


Go 


= n-i/3(2u;)-i/3 , 


E = 


Ec,o + 


n4/3A(2u;)V3 
16^ ^ 


,(10-20) 


iv) V = 


w <0; 


Go 


= 7,-l/3(_^)-l/3 ^ 


E = 


Ec,0 - 




(10-21) 



Here, the quantity Sc,o is defined by Ecfl = Ti-ciqo^Po)- The quantity p {= yo = Vo) 
is obtained from Go with the use of the relation Go = (V(l/2) + p^ + pf- 

For various potentials, let us investigate the energy expectation values under the 
above initial conditions : 

(1) The harmonic oscillator potential uP'Q^ 12 : In this case, v = co'^ and w = are 
obtained. We then have p = {l/2^/2){l/^/Lo — \/oj). Thus, the total energy E is 
expressed as 

E=^pI + + hiu^ . (10-22) 

(2) The Morse potential Tyo[exp(-2//Q) - 2exp(-/iQ)] with Wo > : In this 
case, V = 2p,'^Wo[2 ex.p{—2pqo) — ex.p{—pqo)]. We will show the results of the order 
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for Gq. If qo < log2/;U (case(i)), the total energy E is expressed as 



E = ^Po + ^o[exp(-2/xgo) - 2exp(-/xgo)] + /i/uy -^[2exp(-2^go) - exp(-/xgo)] 2 

(10-23) 

li qo = Pq = 0, we have the minimal classical energy 



Eo = -Wo + nf,^^ . (10-24) 

This energy should be compared with the exact eigenvalue of this quantum Hamil- 
tonian which has the form 



_gexact = 



V 2 



(n = 0,l,2,---) (10-25) 



If n = 0, in the order of , the eigenvalue ( 10-25| ) is reduced to the result ( 10-24 ). 
This means that, even in the limit of ^ — > 0, our formalism automatically contains 
the zero-point energy of order of . 

(3) The upside-down harmonic oscillator potential : In this case, v = —uP' 
and w = (case(ii)). Then, Gq is equal to l/(2u;). 

(4) The double-well potential XiCf' - jl^ with A > : In this case, the cases 
(i), (ii) and (iii) appear depending on the initial value go- If 90 > 1^1 or < —iij^fZ 
(case(i)), Gq and E are given as 



Go 



2A(3(?2 - ^2) 4A(3g2 _ ^ 



2^2 



E = ]g,U VM + A(3,„^ _ + . . (10.20) 

If — ^/\/3 < qo < /i/\/3 (case(ii)). Go is as follows : 



- V 2A(;.2'_ 3,2) + ^ • 8A(/.2 - 3,2)2 • (10-27) 
If qo = ±/i/\/3. Go and E are given as 

Go = n-i(2A)-^ , E = ^pI + V{qo) + nl^{2X)^ . (10-28) 



We can conclude that the conditions in Eq. (10-14) are plausible for describing the 
classical motion in quantal system in our squeezed state approach. The time- 
evolution of the dynamical variables in the squeezed state, which governs the time- 
evolution of the quantum state, is presented in Ref.^2|) in the case of V{Q) = 0, 
V{Q) = 7Q and the upside-down harmonic oscillator potential. 

Before closing this section, we give two comments : First, our squeezed state 
formalism for quantal systems automatically contains the quantum effects beyond 
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the order of h. If we take the limit of ?i — > 0, the WKB result is reproduced. 
Further, we can show that the Maslov phase which appears in the Bohr-Sommerfeld 
quantization rule has the geometrical property such as the Berry phase in the limit 
of 7i ^ in our formalism. Of course, we can take into account of the higher order 



contributions of h. The details can be found in Ref.36). Secondly, we point out that 
our squeezed state formalism can be easily extended to the case of the scalar field 
theory. The details can be found in Ref.33). 



§11. Coherent and squeezed states in the su(2)-boson model 

We extend the squeezed state approach developed in §10 to the su(2)-algebraic 
model. The sti(2)-algebraic models such as the Lipkin model and the pairing model 
are regarded as the schematic models which describe many-nucleon and/or fermion 
system. The S'u(2)-algebra can be expressed in terrns of two kinds of boson operators, 
which is known as the Schwinger boson realization. In this sense, we call the su{2)- 
model treated here the sn(2)-boson model. Focussing on the least quantal effect, in 
this section, we eive an analysis of the sti(2)-boson model based on the coherent and 
squeezed state. E3) 

We introduce two kinds of boson operators a and h. The generators of the 
stt(2)-algebra are expressed by these boson operators, which have already given in 
Eq. ( p-ll ). The Casimir operator -r5u(2) has also been defined in Eq. (2^) with (2-5). 
The magnitude of the su(2)-spin, S", as a mean value is given by 

{w\S\w) = S, (11-1) 

where \w) is a certain state. If \w) is an eigenstate for 5", then the expectation value 
of the Casimir operator is obtained as {w\rg^(2)\'w) = -S*^ + fiS. The second term 
may be regarded as the quantal effect because the mean value of S is 5. In general, 
this relation is not satisfied. 

First, let us adopt a coherent state as \w). The investigation is expected to be a 
powerful help to the classical description of the sii(2)-boson model. As a preparation, 
we introduce the new sets of boson operators : 

c=Ua-Vb, d = V*a + Ub. (11-2) 

Here, U is real and V is complex and they obey the condition 

U'^ + V*V = l. (11-3) 

Then, (c, c*) and {d,d*) are independent boson annihilation and creation operators, 
respectively. Further, we introduce the boson operators (C, C*) and {D,D*) which 
are defined as 

C = c, D = d-exp(-ie/2)xl^ . (11-4) 

V n 

The boson operator D expresses the displacement from the complex value e^"^^ 1"^ ^la jfi. 
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The coherent state |c) has the following standard form : 

|c) = exp(e-^'^/2 • {aa* + /36*))|0) , 

iVe = exp(-(|a|2+/32)/2) . (ll-S) 

Here, a is complex and /3 is real with the help of the introduction of the phase angle 
X- Let a, P and x correspond to 

a = \l^-V, p = ^.U, x = 0. (11-6) 

Here, V, U and 9 have been introduced in Eqs. ( |11-2| ) and ( |11-4| ). Then, the coherent 
state |c) is rewritten as 



Nc = ex.p{-(7/h) . (11-7) 

The state |c) is also the coherent state of c and d with the eigenvalues and 
exp(— • ^2a/h, respectively. Therefore, |c) is the vacuum for C and D : 

c\c)=0, d|c) = exp(-i(9/2) • y— ; C\c) = D\c) = . (11-8) 

It is noted that the coherent state |c) is specified by four variables 6, a, V and V* . 
Let us investigate the quantal fluctuations in the coherent state |c). The squares of 
standard deviations for the quasi-spin operators are calculated as follows : 

{ASD, = {ASl\ = {ASl\ = ^ , (11-9) 

where {ASl)c = {c\Sl\c) - {c\S^\c)'^ etc., and 5^ = (^+ + S_)/2, Sy = (5+-5_)/(2i) 
and Sz = Sq. It should be noted that the standard deviations are fixed in the coherent 
state. Further, it is interesting to investigate the quantal fluctuations by means of 
the expectation value of the Casimir operator rg,^(^2)- We can calculate it as 

(c|An(2)|c) = {clS"^ + nS\c) 

= S^ + hS + ^-{l-W,) , (11-10) 
where Wc is related to the standard deviation for S, i.e., {AS'^)c = {c\S'^\c) — {c\S\c)'^ 

"'--T^ (u.u) 

Here, for the coherent state |c), we obtain Wc = 0. In other words, {AS'^)c is equal 
to hS/2. Namely, it is of the same order as that of HS, which is the exact quantal 
fluctuation. From the above fact, we have to get the conclusion that the coherent 
state cannot give a correct quantal fluctuation with respect to the Casimir operator 
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^su{2)- The fictitious condition Wc = 1 gives us the exact relation for the Casimir 
operator. Therefore, it is inevitable to look for the state \w) which makes the order 
of {AS'^)w less than that of hS. 

In order to overcome the above-mentioned shortcoming of the coherent state, 
we extend the coherent state |c) given in Eq.(ll-5) to the squeezed state \s) in the 
following form : 



\s) = Ns exp(nC'*2 + vD*^) ■ exp e" 



776 



VW+WO ' 

= [(1 + 2^0(1 + 2??* 



v/2(l + 2r]*r]) 



a 



-1/4 . 



expl-- 



(11-12) 

(11-13) 
(11-14) 



It should be noted that the term related to the quadratic for the bosons, which is 
characteristic to the squeezed state, is added. The state \s) contains new complex 
variables and {r]*,r]). Under the condition ( 11-1 ) for \w) = \s), a can be 

related to S in the following form : 



S = {s\S\s) = a + h{CC + r]*r]) 



(11-15) 



First, let us investigate the square of the standard deviation {AS'^)s- The quantity 
{AS'^)s can be calculated as 

{AS^)s = {s\S^\s) - {s\S\sf 



r] T] 



V2 



{r]* +r])r]*r]^/l + 2r]*r]] . (11-16) 



The differential of {AS'^)s for rj and t]* leads us to the following relations for the 
minimum value of {AS'^)s, while we cannot make {AS'^)s vanish any choice of rj and 
rj* for given S and : 



where <^ satisfies the following equation : 



C -2 + 



S 



(11-17) 



(11-18) 



Here, Eq.( |ll-18| ) has a real solution if S/h + (1 - 4^0/4 > 1/2^. Substituting 
the solutions ( |11-17| ) into the relation (11-16), the deviation {AS'^)s is given as 



{AS^ 

ns 



8 V s 



1 + 



8 .Un 



+rc + 2(ro + 



1 n 
6 ' s 



(1 - 4r + 



(11-19) 
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where S and arc regarded as the same order. Thus, up to the order of H^, {AS'^)s 
can be neglected : {AS''^)s = 0. Prom this result, we have 

{s\t^(2)\s) = S^ + hS . (11-20) 

This is the desired result. As was discussed previously, the condition Wg = 1 is 
realized, while Wc = in the coherent state. Thus, the squeezed state \s) gives us an 
appropriate quantal fluctuation for the Casimir operator in contrast to the coherent 
state |c). 

Before investigating the deviations {AS'^)s etc., let us parametrize the state \s) 
in terms of the canonical variables. We have new variables ^* and ^ adding to the four 
variables characterizing the coherent state |c). We impose the following canonicity 
conditions : 



{s\ihdg\s) = S , {s\inds\s) = , (11-21) 
{s\dx\s) = ^, {s\dx*\s) = , (11-22) 

ms) = ^, {s\d,*\s) = ~, (11-23) 



where S is nothing but the expectation value of S and 6 corresponds to its conjugate 
angle variable. A possible solutions of Eqs. (11-22) and (11-23) are obtained as 

T/- / ^ Y TT- ^X*^ ~ * 

V - '^nx*x) ' V 2(5 - 2hx*x) ' ^ " ^ ' ^ - ^ • 

(11-24) 

Then, we can express the expectation values of the generators of the su(2)-algebra 
in terms of these canonical variables : 

(5^), = {s\S^\s) = Qy/S-{P^ + Q'')/4-2nx*x , 

iSy)s = {s\Sy\s) = -P^S - (P2 + Q2)/4 _ 2nx*x , 

{S,)s = {s\S,\s) = -S+ (P2 + q2)/2 + 2nx*x , (11-25) 

where Q = y/h/2-{X* + X) and P = i^/h/2 ■ {X* - X). The expression (11-25) con- 
tains the terms related to the quantal fluctuations, i.e., hx*x. The sti(2)-generators 
are expressed in terms of the quadratic forms for the boson operators. In contrast 
to the case of the coherent state, the forms are not of the normal ordered product 
for the squeezed state \s). Prom this reason, the expressions contain the terms of 
the quantal fluctuations. 

Now we are ready to give the deviations {AS^)s etc., which are expressed up to 
the order of in terms of the canonical variables : 

{ASDs = ^ - [4^2 - Q\4S - (P2 + Q^))w,] 
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+ ^-[AS^-Q\^S-{P'' + Q^))]-x*x , 



{ASl)s = — ■ [452 - P\AS - (P2 + q2))w. 



h 

8S 



■ QP{2S - . -x)^^ + x*x 
+^-[AS^-P\AS-{P^ + Q^))]-x*x , 
{AS% = ^ ■ [452 - (25 - + Q^)fWs] 



~ ■ QP{4S - (P2 + q2)) . _ 1 + 



Here, VF^ is defined by 



Ws = l- = 1 , (11-27) 



where the last equality is valid for the approximation up to the order of as was seen 
in Eq.( ll-l"9| ). In contrast to the case of the coherent state, the above uncertainties 
do not have fixed value hS/2 and Ws is equal to unity even if the variables describing 
the quantal effects, x* and x, are equal to zero. This fact comes from the inclusion 
of the parameters rj and ij* . 

As an example, we sketch our basic idea of a possible application to an su{2)- 
model, that is, the Lipkin model whose Hamiltonian is given as 

H = 2eS, - G{Sl - S^) . (11-28) 

Here, e and G denote positive constant. Main interest is to investigate the time- 
evolution of this quantal system with the least quantal effects. The expectation 
value of this Hamiltonian with respect to the squeezed state \s) is expressed as 

{s\H\s) + 2eS = H, = H,^,i + F,,qi , (11-29) 

^.,ci = e{p^ + g2) + GS{p^ -Q')-j- (P' - Q') ' (11-30) 



Aex*x -^-{P'- - (P2 + q2)) 



~ • [8^2 - 45(P2 + Q2) + (p2 _ g2)2](^* + 3^)^ 1 + 
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G_ 
'2S 



{P^ - Q^)[8S - {P^ + Q^)]x*x 



(11-31) 



The part Hg^d {Hs,qi) denotes the classical (quantal) part of Hg. It is seen from 
(11-29) ^ (11-31) that the c-number Hamiltonian is expressed in terms of the canon- 
ical variables {6, S), {Q,P) and {x*,x). Then, the Hamilton equations of motion are 
written as 



9 = dsHs , S = -deHs = , 
Q = dpHs , P = -dgHs , 



iUx = dx*Hs , 



ihx* 



—dxHg . 



(11-32) 
(11-33) 
(11-34) 



The second equation of Eq. (11-32) shows us that 5 is a constant of motion. By solving 
Eqs. (11-33) and (11-34) with appropriate initial conditions, the time-dependence of 
Q, P, X and a;* can be obtained. Here, let us analyze the Lipkin model under the 
adiabatic approximation like the adiabatic time-dependent Hartree-Fock (ATDHF) 
theory. In the ATDHF theory, the classical Hamiltonian is approximately expressed 
in terms of the quadratic form for the momentum. Following the same viewpoint as 
that of the ATDHF theory, Hg^ci can be expressed in terms of the variables Q and 
P as follows : 



M ^ 



F(Q) = (e-GS)-g2 + ^.Q4 



2{€ + GS) ' 

Further, we approximate i^s,qi in the frame of the term with P° : 



(11-35) 



(11-36) 



In the same way as that in Hg^qi, we have the following approximate uncertainty 
relations : 



4 



(AS'y), - {AS',)s - ^ - {Sx)ld 



Q\^S-Q'){i{x'' -x)X\-^x*x 



1 



X -|- \/ - -|- x*x 



X -\- \\ x*x 
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{2S - QY[i{x* - x)f (^l + x*x^ . 

(11-37) 



Here, (S,),,d = QV S - [P'' + Q'')/i, {Sy)s,ci = -P^^ - (P^ + Q^)/4 and (5,),,ei = 
—S + (P^ + Q^)/2. It is noted that the variable Q should run in the region < 45". 
Our task is to give the initial conditions in order to solve the Hamilton equations of 
motion (11-33) and (11-34). As for the initial conditions for Q and P, they are taken 
arbitrary. It is necessary to give the initial conditions for x* and x. The basic idea 
has been already presented in the previous section, that is, the least quantal effects 
should be realized at the initial time. The idea for investigating the least quantal 
effects in the classical motion of quantized system is given in Ref. ^). Translating 
the basic idea into the present case, the following criteria to give the initial conditions 
should be taken into account for the squeezed state : 

l-ff^qll : minimal , (11-38) 
any two of (D,)° , (Dyf, and (D.f, = , (11-39) 

where the superscript denotes the initial value for each quantity. As can be seen 
in Eq.(ll-37), the criterion (11-39) is satisfied for the case {Dx)^ = {Dz)^ = which 
gives us xq = Xq, that is, xq is real. For convenience, let us introduce the quantities 
g, z and y instead of GS/e, Qq and xq : 

Then, the quantities g and z run from to oo and from to 1, respectively. From 
the criterion (11-38), the initial condition is obtained together with the initial value 
of the energy of the quantal part : 

(1) in the regions < z < zi and Z2 < z < 1, 

g - 2z + 1 - y/g^ - jSz^ - - 4z + 1) +ff - AzJF^ 

y= -5-(4.2-6.+i)+i ' ^''-^'^ 

K^i = , (11-42) 

(2) in the region zi < z < Z2, 



y=, + " + ; . (11.43) 

V -5 - (4z2 _ 6z + 1) + 1 ' ^ ' 



<qi = H^ia • i2z + 1) + l)][-g • (4z2 - 6z + 1) + 1] - (5 - 2z + 1)] . (11-44) 
Here, zi and Z2 are given as 



zi = ^-(g-l + 2sj7g^ + 4^ + 1 cos{4>/^ + 7r/3) 



Z2 = — -[g-l + 2\/7g^ + 4g + l cos(</./3 - 7r/3) 
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7g3 + + 12ff + 2 
cos</> = — , , r, ■ (11-45) 

2( W + 45 + 1)3 

The details are found in Refj35|). It may be interesting to calculate {Dy)^ which is 



given by (Dy)^ = S'^ • f{z) where f{z) = 4z(l — 2z)'^{\ — z)y'^. It can be shown 
that the factor f{z) is rather small with some exceptions. The minimal uncertainty 
for the quantity (Dy)^ is satisfied at rather high accuracy. The detailed numerical 



estimation is given in Ref.35). 



§12. Quasi-spin squeezed state approach to many-fermion systems 



A possible classical description of many-fermion systems such as nucleus is given 
by the TDHF theory E3^. In this theory, the Slater determinant is usually adopted 
as a trial state for the variational procedure. The Slater determinant is essentially 
a coherent state which yields a minimal uncertainty relation where each uncertainty 
is fixed. Therefore, in order to take into account of appropriate quantal effects in 
many-fermion system in terms of classical mechanics, we need to proceed beyond the 
usual coherent state approach. 

In our formalism developed in the preceding two sections, the boson squeezed 
state, which is introduced as a trial state in the time-dependent variational approach 
in the quantum mechanical or many-boson systems, has been originally constructed 
as an extension of the familiar boson coherent state. In many-fermion systems, our 
first task is to define an appropriate trial state similar to the boson squeezed state. 

We formulate our idea for the sn(2)-algebraic models, namely, a single-level 
model with pairing interaction and a two-level model with particle-hole interaction. 
In these models the trial state of the time-dependent Hartree-Fock or Hartree-Fock- 
Bogoliubov theory is reduced to the su(2)-coherent state. With the aim of going 
beyond this theory, we consider a possible extension of the sti(2)-coherent state. 
This new state, which we call "a quasi-spin squeezed state," is used as a trial state 
in the variational method, and we show in this section that the quasi-spin squeezed 
state approach is workable and successful. 

First, we investigate a simple system: N identical nucleons are moving in a 
spherical orbit with the pairing interaction. The single-particle state is specified by 
a set of quantum numbers {j,m), where j and m denote the magnitude of angular 
momentum of the single-particle state (the half-integer) and its projection to the z- 
axis, respectively. The annihilation and the creation operators for the nucleon with 
m, Cm and cj^, satisfy the usual anti-commutation relations : {cm,c'^'} = ^mm' and 
the other combinations are equal to 0. The Hamiltonian is written as 

H= jZ hec*^Cm-\h^G ^ (-y~"^c*^c*_^ i-y~"''c-m'Cm' 

m=—j m=—j m'=—j 

= 2e{So + Sj) - GS+S^ , (12-1) 



where he and h^G denote the single-particle energy and the force strength, respec- 
tively. Also, Sj = hQ/2 where Q is the half of the degeneracy : Q = j + 1/2. This 
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model has been used in §4 as an example of the su(2)-algebraic models. Here, we 
have defined the quasi-spin operators (50,5-1-) as 

% . - h ■ 

m m 

So = ^{T.^^^rr. -n) = |(iV - Q) , (12.2) 

m 

where is a number operator defined by = c^Cm- The set (12-2) is governed 
by the su(2)-algebra : 

SJ\ = 2hSo , [So, S±] = ±hS± . (12-3) 

As is well-known, this model can be solved exactly. The exact eigenvalue Hqi in 
the ground state for this Hamiltonian is given by 

Hqi = {q\H\q) = 2eA - GA(2Sj -A)- TiGA 

= H,i - hGA (= Sqi) . (12-4) 

Here, A = hN/2 in which A^ is the particle number. The state \q) is expressed in the 
following form : 



1 



JV 



Let us investigate the state \w) which gives a classical counterpart for the quasi- 
spin. The state \w) should satisfy the exact relation for the Casimir operator rgy^^2)- 
We set up the following conditions : 

= {w\fg^(^2)\w) = {w\S^ + ^{S-S+ + S+S-)\w) 

= S] + hSj , (12-6) 

Nu, = {w\N\w) = N . (12-7) 

First, we consider a quasi-spin coherent state |c), which is generalized from the boson 
coherent state : 

|c) = C-V2.expf^.%'| 10), 



u h 



k 



(U-^)^ , (12- 



where V is complex and U is real parameter and the relation C/^ + = 1 is 
satisfied. This state is known as the su(2)-coherent state which corresponds to the 
BCS ground-state. It has been regarded as an approximate state for \q) in Eq.( |12-5D . 
This state is a vacuuni state for the Bogoliubov-transformed. fermion operator ' 

am = Ucm - V{-y-^c*_^ , am\c) = . (12-9) 
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By the use of this operator, let us define new operators : 

"(12- 10) 

These operators satisfy the fohowing commutation relations : 



[B, B*] = l- -|-M , [M, B*] = 2B* , [M, B] = -2B . (12-11) 
2Sj 

It should be noted that the state |c) gives the following expectation value for the 
quasi-particle number operator M : n = {c\M\c) = 0. The expectation value of the 
Hamiltonian is expressed as 

i^ch = (clinic) = i^cl - no A + Fch , 

F,^ = —A{2S^-A). (12-12) 

The expectation value is clearly different from the exact energy eigenvalue. Fur- 
ther, it is interesting to investigate the uncertainty relations. Let us introduce an- 
other stt(2)-generators, (s±,so), which consist of the quasi-particle operators a-m 
and instead of Cm and in Eq.(12-2) : s+ = (V^) Z^m(~V~"^^mS-m' ^- = 
ih/2)J2mi-y~'^0'-mam and sq = {h/2) - (Em«mam - ^2)- Then, the following 
squares of the standard deviations for Sx,y,z-, which are defined as Sx = (s+ + s_)/2, 
Sy = (s+ — S-)/2i and Sz = sq, are obtained : 

{AslU = (/^4)ch = ^ , {AslU = , (12-13) 

where {As^)ch = (c|s^|c) — (c|sa;|c)^, etc., and (cls^lc) = —Sj. The above relations 
mean that the state |c) is governed by the minimal uncertainty for the quasi-spin 
operators s^, Sy and Sz, namely, (Z\s^)ch(^Sy)ch = {^"^ /^){c\sz\c)^, etc. Also, each 
standard deviation has a fixed value. We see that the above-mentioned situation is 
similar to the case of the boson coherent state. 

In order to relax the expectation value, n, of the quasi-particle number operator, 
M, and the fixed standard deviations for the quasi-spin operators, we here introduce 
a new state j,s), which we call the quasi-spin squeezed state, as an extension of the 
quasi-spin coherent state |c). The new state is defined as 



|,) = C.-'/^.exp(l.^.B-) 



fc=0 p=0 

For the state \s), we have 



C 

n = {s\M\s) = 4z-^ , (12-15) 



C. 



s 
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where the prime denotes the first derivative with respect to z. The dynamical vari- 
ables in the quasi-spin squeezed state are V, V* , v/u and v* /u. For the canonicity 
conditions for the present variables, the following relations are adopted : 

{s\ihd^\s) = A , {s\ihdA\s) = , 

{s\ind^\s) = hX , {s\indx\s) = . (12-16) 

We then get possible solutions for the above canonicity conditions : 



V = VX- exp(-2z(^ + (j))) , 



The explicit form of A is found in Ref. |37| ). Here, it is easily shown that A is related 
to the quasi-particle number : 

1 71 

X = -{s\M\s) = - . (12-19) 

If A = 0, the quasi-spin squeezed state is reduced to the quasi-spin coherent state. 
Thus, we see that the variables {(j), A) describe the quantal fluctuation and the vari- 
ables A) mainly describe the classical motion. 

The expression of the factor (A-|- 1)(1 — A/i7)(l — (A-|- 1)/17) in is interesting. 
If f] is even or odd, then the maximum of A should be 17 or 17 — 1, respectively. This 
realizes the Pauli principle. On the other hand, if A <C 1?, is approximated as 
(A-l-1). In this case, we have tt^ — |up = 1, which leads us to the boson approximation 
for B and B* . 

Hereafter, the discussion is restricted to the first order quantal fluctuation. In 
this approximation, it is enough to adopt the following form for V, U, v and u : 



V = y^e-"^ , U = ^l-^^, v = ^/Ae-2^(^+^) , u = VT+X . (12-20) 

Then, the expectation value of the Hamiltonian with respect to the quasi-spin 
squeezed state is approximately given as 

i/sq = {s\H\s) = H,i - fiGA + , (12-21) 



/(A, 0) = 1 + 2A + 2y^ A(A + 1) - cos 20 . (12-22) 

Here, Fch has been defined in Eq.(12-12). In this approximation, the squares of the 
standard deviations for the quasi-spin operators Sx^y,z are given as follows : 





)sq 


_ T^Sj 
2 




)sq 


_hS, 
2 




)sq 


= , 



2sJX{X + l) -cos[2(^' + , 

(12-23) 
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where (Z\s^)sq = (s|s^|s) — {s\sx\s)'^, etc., and (sisals) = —Sj. 

Our problem is to solve the equations of motion derived from the time-dependent 
variational principle : 

• ^ _dHsq ■ ^ dHsq 

d0 ' OA ' 

« = ''^ = (1^-^^) 

The initial conditions for A and (j) which describe the classical image of quantal 
fluctuation are obtained under the similar criterion to that given in the case of 
quantum mechanical system in §10 or the SM(2)-boson model in §11. Namely, at 
the initial time, the minimal uncertainty should be realized. Further, the energy 
originated from quantum fluctuations is minimal at the initial time. 

• = Q)' • \{s\s,\sf\^ , fiXoAo) ■■ minimal , (12-25) 

where the superscript and subscript "0" denotes the initial value of corresponding 
variable. These criteria lead us to the following initial conditions : 

</'o = ^-^o, (A; = 0,±l,±2,---) , x^ = U \ -l) . (cos2,^o<0) 
2 2 V I sm 20Q I / 

(12-26) 

If cos 2^0 > 0) then Aq = 0. For the solution ( 12-26|) , i/gq is given in the form 



hG 

= Hd - fiGA + — • A{2Sj -A)-\ sin 2^q\ . (12-27) 
2b j 

It may be interesting to see that if |sin2^o| vanishes, Fsq = and if |sin2^o| is 
equal to 1, Fgq = hG/2Sj ■ A{2Sj — A). The former and the latter correspond to the 
situations given by \q) and |c), respectively. Therefore, by the parameter |sin2^o|i 
we can reproduce the intermediate situation. Concerning the minimal uncertainty, 
we have 



= Vl + |sin2^^j±Vl_-^in2go| _ 
v/2|sin2^o| 

The relation (12-28) shows us that at the initial time the magnitudes of (Z\s^)gq and 
{Asy)gq are controlled by the quantity sm2^Q, which is the same situation as that in 
the case of Fgq- Under the above initial conditions, we obtain the following solutions 
for the Hamilton equations of motion : 

A{t) = A , 

^(t) = 2uJsqt + ^0 , 

A(t) = l^j A\2Sj - Af\sm2'Po\t^ 
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+^A{2Sj - A)\ cos 2^o\t + ^ I 



1 



sin 2#o| 



-'J 

= i • cos-i x(i) , (12-29) 
where Wgq and x(f) are defined as 

= e - G{S, - ^) - ^ • (1 - | sin 2^o|) - ^| sin 2^o| , 

^ |sin2^o|-(l + 2A(t)) ' 
"^^^ 2^/A(0(A(i) + l) ^ ^ 

Next, let us investigate another many-fermion system which consists of two- 
energy levels with same degeneracy 2f2 (= 2j + 1). The two-body interaction is 
active only for a particle-hole pair with coupled momentum being 0. This model is 
called the Lipkin model whose Hamiltonian is written as 

H = 2eSo-^G{Sl + S^_) , (12-31) 

where 2fie is the level spacing between lower and upper energy levels labeled hy g = 1 
and 2, respectively. We have already used this model in §§4 and 11 as an sii(2)-boson 
model. However, in this section, we deal with this model as an example of many- 
fermion system. Thus, the quasi-spin operators consist of the fermion operators, 
which are defined as 

. . - ^ - n ^ 

S+ = h ^ C2^Ci^ , = h ^ cl^C2m ) So = — ^ (C2m'^2m — Cim^lm) , 

m=—j m=—j m=—j 

(12-32) 

where c*^ and Cgm are fermion creation and annihilation operators. Here, the 
above quasi-spin operators satisfy the commutation relations of the sit(2)-algebra 
: [5"+ , S-] = 2hSo and [Sq , S±] = ±hS±. Hereafter, we concentrate on the case 
N = 2Q where N represents particle number. 

The stt(2)-coherent state is adopted as a trial state in the time-dependent Hartree- 
Fock theory. This state is given as 

1 / c \ 

(12-33) 

with C2m|0ph) = Ci^lOph) = 0. This coherent state is a vacuum state with respect to 
the unitary-transformed fermion operators as 

. _ 1 . Z ^ 

"2m - (1 + |_^|2)l/2'^2m - (■]_ _^ |_^|2^1/2'^1'" ' 

Z* 1 

^1"^ = (l + |^|2)l/2^2m + (1 + |_^|2)1/2^1"^ (12-34) 

with a2m|c) = Oimk) = 0- 
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Our aim is to extend the coherent state to the squeezed state in order to go 
beyond the usual Hartree-Fock (HF) approximation. Now, let us introduce a quasi- 
spin squeezed state from the sn(2)-coherent state in a manner similar to that for the 
case of the pairing model discussed previously. The new operators which consist of 
the unitary-transformed quasi- fermion operators are introduced like ( 12-10| ) : 



m=—j ' m=—j 

j 

M = ^ (a2„a2m + aimOlm) • (12-35) 
m=-j 

The operator M is nothing but the number operator of the quasi-particles and holes. 
The commutation relations among these operators have the same forms as Eq.( 12T"l| ) 
except for the replacement of Sj to f2. The quasi-spin squeezed state in this model 
is defined as 



• exp I B ^ 

2 u 



Let us investigate the classical counterpart given by the quasi-spin squeezed state. 
We impose the canonicity conditions so as to transform the variables, which are 
contained in the |s), from Z, Z* , v/u and v* /u to b, b* , fi and (3* : 

/l^l\ "^h* I \ ^ \ \ /l^l\ /l^l\ ^ R 

^%"^^ = 2^ ' ^'\W>^ = -2^-^ ^^"9^"^^ = 2^ ' ^^"9^"^> = -2^- 

(12-37) 

Under these conditions, the sets of variables (6, b*) and (/?, /?*) are those of canonical 
variables, respectively. Solving the above canonicity conditions, we can express the 
expectation values of the quasi-spin operators in terms of 6, 6*, (3 and (3* : 



{s\S+\s) = Vnb*^2n - h\b\'^ - 4n|/3|2 = , {s\So\s) = -Q+h\b\^+2n\l3\'^ . 

(12-38) 

These expressions remind us of the Holstein-Primakoff boson representation for the 
SM(2)-algebra. Further, we can calculate the energy expectation value with respect to 
the quasi-spin squeezed state. The result is almost identical with the exact eigenvalue 
for the ground state in the wide region of the coupling strength G. The details are 



found numerically in Refs. 35) and|3^). 

Finally, let us investigate the implication to the RPA.eHP We now set the operator 
as a pure boson operator B. Further, under this replacement, the number operator 
of the quasi-particles and quasi- holes, M, should be replaced so as to satisfy the 
original commutation relation ( |12-llD . Namely, 

B — >B , B* — >I3* , [B , I3*] = l ; M — >M = 2B*B . (12-39) 
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Also, the quasi-spin squeezed state \s) has to be replaced into the usual boson 
squeezed state |6) : 



|6) ^ r-^/^xp ( - • 



1 V 



*2 



) , (12-40) 

\ ^ a, / 

where the relation ti^ — |up = 1 is satisfied. As is well known, the boson squeezed 
state is a vacuum with respect to the unitary-transformed boson operator : 

6|&)=0, b = uB-vB* , [b,b*] = l. (12-41) 

Let us consider to diagonalize the model Hamiltonian with respect to the boson 
operators b and 6*. The Hamiltonian ( 12-31 ) is expressed in terms of B* and M 
as 

H = E,h + EuM + E2oB* + E;oB + EMMM^ + E22B*B + {E^oB*^ + EmBM + h.c.) , 

(12-42) 

where Ech = (c|-ff|c), Eij denote c-number coefficients and h.c. means the Hermitian 
conjugate. In the HF approximation, £"20 = E2Q = are demanded because B* is a 
particle-hole pair-creation operator. Now, let us take into account of the higher order 
terms, which are neglected in the HF approximation. Following the replacement 
( 12-39| ), neglecting the terms of third and fourth order with respect to the boson 
operators B and B*, and transforming i3 to 6 in ( |12-41 ), we obtain the Hamiltonian 
as 



B 

Hf, ~ ^ch + a\v\'^ + !^{uv + uv*) + 



+ 



a(n^ + |^;|^) + P{uv + uv*) 



auv H — [u + V ^ 



b*b 



b*^ + h.c. 



(12-43) 



Here, if we define x = fi^Gf2{l - l/2f2)/{he), then a and /? are given as a = 2he for 
X < 1 or 7ie(3x^ — l)/x for x > 1 and (3 = —2hex for x < 1 or —he{x'^ + l)/x for 
X > 1, respectively. We demand that the terms which is proportional to 6*^ and 
should be zero. As a result, in the limit of 2i7 » 1, u and v are determined as 



1 



1 + 



a 



The Hamiltonian is finally expressed as 
Hb ^ E^pA + hujfb , 



^RPA = ^ch + ^(\/a'-/3' 



a] 



1 



a 



- /32 



(12-44) 



a2 - /32 . (12-45) 



Here, -Erpa is identical to the ground state energy of the RPA calculation. Namely, 
the ground state correlation is included automatically in our approach. The boson 
operator b corresponds to the phonon operator. Thus, our quasi-spin squeezed state 
approach includes the RPA in a certain limit. 
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§13. Concluding remarks 

On the basis of tlie celebrated boson expansion theory^ proposed by Marumori, 
togetlier with Yamamura (one of the present authors) and Tokunaga, we have devel- 
oped the idea of the boson mapping method, which is contained in the Marumori- 
Yamamura- Tokunaga boson expansion theory, to describe the thermal states in quan- 
tum many-body systems. The original boson mapping method in which the fermion 
space is mapped to the physical boson-space is not restricted to the many-fermion 
systems such as nucleus. We have used this method in §§2 and 3 to give the classical 
counterparts of the su{2)- and the sn(l, l)-algebras realized by the Schwinger boson 
representation. As a result, it has been shown that the classical counterparts of 
these algebras realized in terms of the Schwinger bosons could be obtained in the 
same forms as those of the Holstein-Primakoff representations of the su{2)- and the 
su{l, l)-algebras, respectively. With the aim of describing the thermal behaviors of 
the quantum many-body systems, we have learned in §§4 and 5 that the thermal 
properties appear with the su{l, l)-algebraic features in the "damped and amplified 
oscillator" model. With this insight, in §6, we have considered the sii(2)-spin system 
and our idea has been applied to this model. The ■SM(2)-spin model is often regarded 
as a schematic model of nucleus. It has been indicated that the sti(2)-spin model 
shows the su{l, l)-like behavior with the help of a certain kind of the coherent state 
giving a classical counterpart. Further, we have treated the su(2)-spin system in- 
teracting with an external harmonic oscillator. With the help of the MYT boson 
mapping method again, we have shown that the ■sm(1, l)-like behavior is realized in 
the physical space spanned by the boson mapping from the original space. Thus, it 
becomes possible to describe the thermal behavior in the system composed of the 
■su(2)-spin interacting with the harmonic oscillator. It is seen in §7 ~ §9 that our 
formalism is workable to describe the thermal behaviors in this model. 

Another powerful method to describe the quantum many-body systems in the 
microscopic way is the time-dependent Hartree-Fock theory. To investigate the small 
amplitude nuclear collective vibrational motion, the idea of the TDHF theory was 
proposed by Marumori. 0* In the TDHF theory, an approximation of the small am- 
plitude fluctuation around a static Hartree-Fock field results in the random phase 
approximation. This approximation is not enough to investigate the large amplitude 
nuclear collective motion. To go beyond the RPA, the TDHF theory devised by the 
canonical form was proposed by Marumori, I12P together with Maskawa, Sakata and 
Kuriyama (one of the present authors), in which the collective coordinate was ex- 
tracted self-consistently from many-nucleon degrees of freedom. The time-dependent 
variational approach with the canonical form developed by Marumori et al. is not 
restricted to treating the many-nucleon or -fermion systems. Also, it is not restricted 
to using the Slater determinant, which is essentially a coherent state, as a trial state 
for the time-dependent variational procedure. We have extended this formulation to 
quantum mechanical systems and quantum many-boson systems in order to describe 
the time-evolution of these systems in terms of the classical mechanics. The key idea 
is to use squeezed states as trial states in the time-dependent variational procedure. 
As a result, it has been shown in §§10 and 11 that the higher order quantal effects 
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than the order of fi^ have been automatically included in our squeezed state formal- 
ism. Further, our idea is applied to the quantum many-fermion systems governed 
by the sti(2)-algebra. We have defined and introduced the quasi-spin squeezed state 
to go beyond the usual TDHF theory in §12. It has been shown that, in the single 
orbit shell model with pairing interaction, the BCS result and the exact one for the 
ground state have been connected by the use of this new state. Also, in the Lipkin 
mode, we have shown that the quasi-spin squeezed state approach could reproduce 
the result obtained in the RPA under a certain limit and, further, seemed to give us 
a possible way to go beyond the RPA. 

The boson expansion theory and the TDHF theory in the canonical form are 
powerful theories to investigate diverse quantum many-body systems in various situ- 
ations. It is expected that these theories will be developed in and successfully applied 
to quantum physics in addition to the microscopic studies of nuclear collective mo- 
tions. For example, we can find such studies in Ref. 43). 
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